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This paper develops the basic sets of equations which lead to the conservation laws describing 
collisionless plasma shock waves. We discuss the evolution of shock waves by wave steepening, derive 
the Rankine-Hugoniot conditions for magnetogasdynamic shocks, discuss various analytical models 
of shock formation, and discuss the basic instabilities which may become important in collisionless 
shock physics. We then present a survey of the theory of anomalous resistivity in the quasilinear limit 
and beyond and discuss mechanisms of shock particle reflection as far as they have been investigated 
qq in the published literature. The content of the chapter is the following: 1. Wave steepening, 

describing simple waves and steepening due to nonlinearity, balnced by dissipation in Burgers' 
shocks, by dispersive effects in the Korteweg-de Vries equation, the Sagdeev-Potential method, 2. 
£SJ Basic equations, presenting kinetic theory and the transition to moment equations in the fluid 

description, 3. Rankine-Hugoniot relations, giving the jump conditions across shocks, explicit MHD 
solution for perpendicular shocks and parallel shocks, and high Mach number conditions, 4. Waves 
and instabilities, giving the general dispersion relation, describing low-/3-shocks, whistler and Alfven 
shocks, the various shock-relevant instabilities, 5. Anomalous transport for the various electrostatic 
■ wave-particle interactions, general description of anomalous resistivity, shock particle reflection from 

I potential and specularly, hole formation, 6. Briefing on numerical simulation techniques, giving a 

short idea on this important field and its methods. 

C| PACS numbers: 

O | Keywords: 

6 

, I. WAVE STEEPENING 

Shocks have a certain width A and a certain jump in density N, temperature T, pressure P and magnetic field B 
across this width from a given upstream value to a downstream value. This jump is by no means infinitesimal. At the 
\^ contrary it is usually several times the upstream value in magnitude. Thus, looked at as a wave, a shock is a highly 
nonlinear wave structure of wavelength A with amplitude that cannot be neglected compared with the upstream 
^SJ value. Therefore, the basic equations describing a shock cannot be linearized as is usually done in considering wave 
jy-^ phenomena. These equation must be solved in their full nonlinearities. This, however, is barely possible and can be 
done analytically only in very rare cases which usually are not of interest. On the other hand, shocks evolve inside 
the plasma from small disturbances. It is thus reasonable to ask for the nonlinear evolution of such a small harmonic 
disturbance in order to learn, how a disturbance can evolve into a very large amplitude shock ramp. 

> 



X 



A. Simple waves: Steepening and breaking 

The simplest way to do this is to consider the evolution of so-called simple waves [see, e.g., [56] . Simple waves are 
one-dimensional sinusoidal disturbances of the plasma velocity of the form V(x) = Asinkx moving on the plasma 
background at speed c. The total derivative of this disturbance is given simply by 

dV dV Tr dV n 

d^^ + ^ = ° (1) 

and in the absence of any forces and friction is assumed zero. Hence this is the equation which describes the evolution 
of the disturbance as long as no friction or other force comes into play which in the initial state is a reasonable 
assumption. 
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FIG. 1: Schematic of the steepening and wave breaking phenomenon, illustrated for three successive times ti, t2, £3. At t2 the wave has 
steepened to maximum, and in t-g, it collapses in the absence of any retarding effects. 
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FIG. 2: Steepening and final breaking of an initial sinusoidal (simple) wave. Left: Initial wave form in wave frame showing the 
nonlinear action of the wave on its own shape. Right: Calculation of wave form steepening in a shallow fluid of depth h [after [ST] . The 
wave profile is shown at the initial time to, intermediate time t\, and breaking time tj, when the wave starts turning over. Steepening of 

the profile is well expressed. 



Let us now investigate how such a disturbance will evolve when it propagates through the plasma. Clearly during 
propagation the main effect on the shape of the disturbance arises from the second nonlinear term which can be 
written as Vkcoskx. Inserting for V this becomes |sin2fca;. Hence, the nonlinear term in the above equation |l| 
generates harmonic sidebands of half wavelength and half the amplitude of the original wave. These waves, by the 
same mechanism, also generate sidebands on their own now at quarter original wavelength and amplitude, and so on 
with increasingly shorter wavelengths. The total amplitude is the superposition of all these sideband harmonics 

V k {x) = Y J Ji^2 l kx, (I = 0,1,2,...) (2) 

1 

All these waves and sidebands propagate at the same velocity c. This can be easily seen when, for propagating waves, 
replacing kx — > (x — ct) in the above expression. They locally superpose and add to the wave amplitude. Because 
ever shorter wavelengths contribute, the wave steepens until the gradients become so short that other processes take 
over. If this does not happen, the wave will turn over and break. 

This is illustrated in Figure[T] and [2] in the co- moving frame of a sinusoidal wave. The left part of the figure shows 
the nonlinear mechanism. Since the velocity is largest at the maxima it speeds up the motion of the maxima with 
respect to the remaining parts of the wave profile. Moreover, the action on the positive and negative maxima are 
oppositely directed, and the wave starts forming a ramp corresponding top a shock front. This happens at time t^. 
For times t > % the wave will turn over and collapse. This can be prevented only by additional processes which set 
on when the wavelength of the ramp becomes so short that in Eq. ([IJ terms of higher gradients in the velocity must 
be taken into account. 

Eq.Q can in fact been understood as the lowest order equation describing the evolution of a wave packet of wave 
number k. In general in the wave frame of reference its right-hand side is a function F(V) that can be expanded 
with respect to V. The first higher order term in this expansion turns out to be second order in the spatial derivative 
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V = xd x , where x is the unit vector in the direction of x. The next higher order term is third order in V, and so on. 
Up to third order the resulting equation then reads 

— + V— - — D— - p— + (3) 
dt dx dx dx dx 3. 

The first term on the right is a diffusive term with diffusion coefficient D(x). The third term with arbitrary coefficient 
(3 is the lowest order contribution of wave dispersion to the evolution of the wave shape and amplitude. This can be 
most easily seen when taking the linearized equation, assuming V to be a small disturbance only of the wave speed 
c and neglecting the nonlinearity in the second term on the left by approximating it with cW, and subsequently 
Fourier analysing for a harmonic perturbation V — Aexpi(kx — u>t) with wave number k and frequency u). This 
procedure yields (for constant D) the following dispersion relation 

oj-kc+k 3 f3 = ~ik 2 D (4) 

On the left of this equation there is the relation between the frequency u) and wave number k, while on the right 
appears an imaginary term that depends on the diffusion coefficient D and square of the wave number k. Imaginary 
terms in frequency imply damping. Hence, as we did express above, the second order spatial derivative term in 
Eq. (J3j> corresponds to diffusive dissipation of flow energy, while (for real (3) the third term in Eq. (|3| causes the wave 
to disperse, i.e. waves of different wave- numbers, respectively different wavelengths, propagate at different phase 
velocities. 



B. Burgers' dissipative shock solution 



Returning to the nonlinear equation ([3| we can thus conclude from this analysis that a wave if not breaking will 
steepen so long until either the diffusive or the dispersive terms on the right start competing with the nonlinearity 
in which case the wave may assume a stationary that is balanced either by diffusion or by dispersion. Diffusion 
implies balance by real dissipation with energy transformed into heat, while dispersion implies that the 'dangerous' 
short wavelength waves which cause the steepening either run the wave out or do not catch up with the wave when 
the slope of the wave profile exceeds a certain steepness. In the dispersive case the wave will either exhibit short 
wavelength fluctuations in front of the steepened profile or behind it, depending on whether the shorter wavelength 
sidebands are retarded or accelerated. However, we can also conclude that dispersion alone should be unable to 
generate a shock since in the simple form discussed here it does not produce irreversible dissipation and hence no 
heating and increase in entropy. For a shock profile to be created some kind of diffusive process will be necessary. 



1. Stationary Burgers equation 

Eq. ([3]) allows to distinguish two extreme cases. The first case is that of purely diffusive compensation of the nonlinear 
steepening. In this case the dispersive term can be neglected, and one obtains the Burgers equation 

dV Tr dV d 2 V 

-m +v dx- = D o^ (5) 

which is a non-linear diffusion (or heat conduction) equation. In contrast to the ordinary linear heat conduction 
equation the Burgers equation possesses stationary solutions due to the above mentioned compensation of diffusive 
spread by nonlinear steepening. These stationary solutions can be found when transforming to a coordinate system 
moving with the wave by introducing the new coordinate y — x — ct. Then Burgers' equation becomes 

d 2 V dV 

We are interested only in solution which are regular at infinity with vanishing derivatives. Introducing the variable 
V' = V — c the first integral is easily obtained. Integrating a second time the solution found is then 

V , fx-ct\ 

The form of this solution is a typical shock ramp which is displayed in Figure|3] The ramp is sitting on the wave 
velocity c. Its width is A = 2D/c. The shock solutions produced by Burgers' equation are thus propagating non- 
oscillatory shocks; they are simple stationary constant amplitude ramps of the sort of tsunamis. We should, however, 
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FIG. 3: The stationary solution of the Burgers equation is a smooth shock ramp of width A = 2D/c, depending on shock velocity c 

and diffusion coefficient D. 



keep in mind that they are produced solely by nonlinear steepening and its compensation through diffusion. When 
the latter is large, the shock will be steep; in the opposite case it will be a flat ramp only, and its relative height is a 
function of its width. 



2. Time- dependence 



Since Burgers' equation is an ordinary diffusion equation it can be solved by the usual methods of treating heat 
conduction evolving from some initial state. Such an investigation is necessary in order to justify that the stationary 
state of the shock ramp Burgers solution can indeed be reached by evolution out of the initial state Vq. To this end one 
transforms the time-dependent Burgers equation (pi through introduction of a new variable <j> via V — —2D[d In 4>/dy] 
into the common form of a diffusion equation d<f>Jot = Dd 2 cf>/dy 2 , which has the commonly known solution 
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The initial disturbance satisfies the condition of convergence f? dy'Vo(y') < const • y for y 
requirement that f_ &y'Vo(y') — < oo as well as the time-asymptotic solution 



(8) 
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The dummy function G(x) is a transformed version of the function <f> that is given by 7r a G(x) = e ®/ 4D f* d?/e v + 

e e/4D j°° (\rje~ 71 . This asymptotic profile Burgers shock ramp solution is shown in Figure 4 The characteristic shape 
of this solution contains a smooth wavclikc increase up to a flat plateau followed by the shock ramp and a smooth 
transition to the undisturbed state. The ramp is moving to the right in the direction of original wave propagation. 
This is seen from the time dependence of the crest of the ramp. Clearly, most shock transitions in space do not 
exhibit this smooth rise inside the downstream region of the shock, indicating that the Burgers solution has pure 
model character which does not really confirm with the plasma reality. 



C. Korteweg-de Vries dispersion effects 



Balancing the nonlinearity with the help of dissipation is one possibility. The remaining possibility which in the 
absence of dissipation becomes the dominant, is balancing nonlinearity with dispersion. In this case we can neglect 
the diffusion term in Eq. ([3| to obtain the so-called Korteweg-de Vries equation |91j 



(10) 



2(Df) 1/2 

Distance y 



FIG. 4: The time-asymptotic Burgers' shock solution Eq. |9| which evolves from the initial disturbance V?0 through steepening and 
dissipativc ramp formation after a given diffusive time at a location which is determined by the time t and diffusion coefficient D. The 
characteristic shape of this solution contains a smooth wavelike increase up to a flat plateau followed by the shock ramp and a smooth 
transition to the undisturbed state. The ramp is moving to the right in the direction of original wave propagation. This is seen from the 
time dependence of the crest of the ramp. Clearly, most shock transitions in space do not exhibit this smooth rise inside the downstream 
region of the shock, indicating that the Burgers solution has pure model character which does not really confirm with the plasma reality. 



Similar to Burgers' equation, the Korteweg-de Vries equation also allows for stationary localized solutions. Such 
solutions are restricted to a finite spatial interval because dispersion does not cause irreversible effects. As before we 
assume that the stationary solution moves a speed c, and we introduce the co-moving coordinate y — x — ct, this 
time measured from the centre of the localized disturbance. The Korteweg-de Vries equation then transforms into 
the third order ordinary differential equation 

,„ ,dV n d 3 V n 

Solution of this stationary equation requires the prescription of boundary conditions at y — > ±00 for which we choose 
V = dV/dy = 0. It can be shown by substitution that the function 



Vk<iv( x — °t) — 3csech 2 



c (x — ct) 
2 



(12) 



This function is a so-called soliton solution; it describes a stationary bell-shaped solitary wave pulse propagating at 
velocity c along x without any change of form. The width of this pulse is A = 1\J j3/c and depends on the velocity 
c and the dispersion parameter (3 in such a way that the faster the pulse moves, the narrower the pulse becomes. In 
addition there is a distinct relation between the amplitude A of the pulse and its width 

A = 2(3(3/A)i (13) 

from where it follows that large amplitude Korteweg-de Vries solitons are fast and narrow. Like in the case of the 
Burgers equation, the Korteweg-de Vries equation is a model equation which results from the dispersive properties of 
nonlinear waves. However, it is interesting that it can be derived for real problems arising in plasma wave propagation, 
and several variants of it have in the past been applied to the plasma. Hence, in describing the nonlinear evolution of 
plasmas it is a more realistic model than Burgers' equation which is simply a consequence of the strongly dispersive 
and practically dissipation-free properties of plasmas. 

The stationary Korteweg-de Vries can have a whole chain of such solitons as solution with the solitons having 
completely different amplitudes and widths. These solitons have the interesting property that they can pass through 
each other during mutual encounters without having any effect on their widths and amplitudes; only the phases 
and spatial positions of the waves from which the solitons form will change during the collision. The question of 
how these chains of solitons are produced is a question that can be answered only when solving the time-dependent 
Korteweg-de Vries problem imposing a certain initial condition similar to that imposed above on the time-dependent 
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FIG. 5: Top: Shape of Korteweg-de Vries solitons for different widths and amplitudes (labelled Ax, A2, A3). Solitons of large amplitude 
are narrower than those with smaller amplitude. Below: In the interaction between two Korteweg-de Vries solitons of different speeds 
only the position changes while the collision has no effect on amplitude or shape. 



Burgers equation. In the case of the time-dependent Korteweg-de Vries equation the solution cannot be found in such a 
simple way, however. Solving it rather constitutes a major mathematical problem which requires solving an equivalent 
Schrodingcr equation [5D]. It turns out that the soliton amplitudes in the chain which solves the Korteweg-de Vries 
equation are related in some manner by an infinite set of invariants of the Korteweg-de Vries equation. In addition, 
the time-dependent Korteweg-de Vries equation also supports wave trails which accompany the solitons forming an 
oscillatory (turbulent) background of spatially dependent amplitudes on which the solitons propagate. 

Clearly, these soliton chains are no shocks; they are wave pulses which after some steepening time evolve into sta- 
tionarity and are completely reversible practically not leaving any effect on the plasma if one neglects the microscopic 
processes which take place in the plasma. But this is precisely the door where the speculation comes in about such 
solitons in collisionlcss plasma being the initial state of the formation of collisionless shocks. Because, if one can 
manage a soliton in the chain to move so fast that its width becomes comparable to the intrinsic plasma scales, then 
the wave field of the soliton should distort the microscopic particle motion causing some kind of dissipation which 
necessarily will turn the soliton into a shock wave by generating entropy and producing a difference between the 
states upstream and downstream of the soliton. The soliton in this case borrows from Burgers shock solutions, and 
mathematically the Korteweg-de Vries equation regains the lost dissipative Burgers second-order term becoming a 
Korteweg-de Vries-Burgers equation. This was, actually, the point |46j made intrinsically in his famous theory of 
shock formation in collisionless plasma. 

D. Sagdeev's Pseudo-Potential 

The Korteweg-de Vries equation is the ideal candidate for introducing one particular notion that has become immensely 
important in soliton and shock research, the so-called Sagdeev potential. The Sagdeev potential is a pseudo-potential 
introduced in order to solve a certain class of nonlinear partial differential equations and to distinguish between solitary 
wave and shock solutions of these equations. This method takes advantage of the similarity of the first integral of the 
particular class of equations to the equation of motion of a hypothetical particle in classical mechanics. Knowledge 
of the Sagdeev potential then reduces the problem of solution to the mere discussion of the behaviour of a particle in 
the pseudo-potential well. 
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V=2c 



FIG. 6: A sketch of the Korteweg-de Vries Sagdeev pseudo-potential. Solutions exist only in the region of S(V) < 0. The maximum 
soliton amplitude is just V = 3c. The minimum potential is at V = 2c. The dotted line shows the path of a pseudo-particle for shock 
formation in presence of dissipation when it "steps down" the potential well to approach the "minimum energy" state. 



The stationary Korteweg-de Vries equation (111 can be directly integrated once. Applying the boundary conditions 
at infinity, the integration constant in the first integration becomes zero yielding the nonlinear second-order differential 
equation 

d 2 V ( 1 \ 

The similarity of this equation to Newton's equation of motion of a pseudo-particle of mass (3 in the force field given on 
the right-hand side is quite obvious. Here V is the pseudo-spatial coordinate, and y is the pseudo-time. This equation 
can be solved by multiplying it with the pseudo-velocity dV/dy, after which the right-hand side can be represented 
as the derivative of the Sagdeev pseudo-potential, S(V), which in this case is a function of the (real) velocity V. It 
becomes explicit in the pseudo-energy conservation law 

after having integrated a second time and again applied the vanishing boundary conditions at infinity. Because the 
left-hand side of this expression is a positive quantity, solutions exist only under the condition that the Sagdeev 
pseudo-potential is attractive, i.e. is negative 

and, as we already know, solutions can exist only in the region of velocity space where V < 3c. Figure|6] shows a sketch 
of the Korteweg-de Vries Sagdeev potential. It vanishes at V = and V = 3c and has its minimum of S = — 3c 3 /2 at 
V = 2c. In terms of energy states a pseudo-particle (soliton) can assume any of the energy levels inside the negative 
portion of the Sagdeev potential. The soliton with maximum amplitude V — 3c is at the "highest" level S — 0. But 
there can be many solitons at this level with amplitudes between and 3c. The "most stable" soliton at the "ground 
state" has minimum Sagdeev potential and amplitude V — 2c, and there is only one soliton with such an amplitude. 
In the absence of dissipation all these solitons are stable. 



The actual solutions ( 12 1 of the Korteweg-de Vries equation can then be found from dV/dy = ^-2S(V)/P by 



simple quadrature, solving the integral 



yo = I dV t (17) 



o 
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We do not further discuss the stationary solution of the Korteweg-de Vries equation; it is but one example which can 
be solved by the Sagdeev potential method. In the literature it has been demonstrated that a very large number of 
other nonlinear problems in plasma related to solitary waves can be treated in the same way, sometimes under much 
more complicated conditions and leading to different types of solitary solutions. 

All these solutions are, however, dissipation free and do not directly lead to shock wave solutions. As [JS] pointed 
out, they will turn into shock solutions whenever anomalous processes at short wavelengths cause the appearance of 
some kind of anomalous dissipation under the ideal conditions of non-collisionality. 

This claim is a most important insight that can, however, be based only on the kinetic theory of the microscopic 
interaction between waves and particles and waves and waves in plasma far from thermal equilibrium, the so-called 
collective processes which dominate the behaviour of high temperature plasmas in which shock waves are at home. 
The dotted arrowed line in Figure[6] shows the presumable "path" of such a dissipative soliton in the Sagdeev potential 
"energy" space. The soliton pseudo-particle will in this case step down the potential, possibly in an oscillatory way. 
In the case when the system is open and energy is continuously supplied it might reach a stationary shock state with 
shock amplitude V s or settle at the ultimate minimum of the Sagdeev potential. The dynamics of this depends on 
the microphysics. 

II. BASIC EQUATIONS 

Before discussing these processes and their relevance for shock wave formation, we need to briefly introduce the 
equations which lie at the fundament of all these processes and to discuss their macroscopic consequences. We 
will in the present chapter distinguish between two approaches to the description of shocks, the theoretical and the 
numerical approaches, respectively. The former deals with the average properties of collective plasma behaviour and 
the investigation of wave growth from an infinitesimal perturbation up to a large amplitude shock, the latter refers 
to the dynamics of macro-particles (as has been described in Chapter 1) and is independent of the average equations 
as it simply solves Newton's equations of motion of the many macro-particles that constitute the plasma in their self- 
consistent fields, where the fields are obtained from Maxwell's equations of electrodynamics, any final shock theory 
must combine both approaches because the fundamental basic equations cannot be solved analytically, while the 
numerical approach provides data which cannot be understood without a follow-up theoretical investigation tailored 
to serve the effects found in the numerical simulation experiments. 

A. Kinetic plasma equations 

Collisionless shock waves represent the final result of collective interactions in which very many particles and in 
addition the full electromagnetic fields are involved. It is thus quite reasonable that they cannot be described by test 
particle theory which considers the motion of non-interacting particles separate from other particles and fields. Test 
particle theory can however be applied if one is not interested in the formation of shocks but instead in its effect on 
small numbers of particles. This is used in the shock theory of charged particle acceleration which will be the subject 
of Chapter 6. 

1. Maxwell-Vlasov equations 

The basic equations on which shock physics is founded are the kinetic equations of a plasma [cf., e.g., HOI HHJ and 
others] or, at the best, some of its simplifications, in addition to the full set of the equations of electrodynamics. 

Since collisions can be neglected, and thus the Boltzmann collision term in the kinetic equations is suppressed, 
these equations reduce to the (non-relativistic) Vlasov-Maxwell set of equations 

ap± „ 

— +vVf ± + -(E + vxB).V/=0 (18) 
at m± 

where F ± (v,x,t) are the electron and ion phase space distributions, distinguished by the respective + and — signs, 
which depend on the six-dimensional phase space velocity, v, and real space, x, coordinates. m + = rm and rn_ = m e 
are the ion and electron masses, respectively; e + — e is the ion charge, e_ = — e electron charge, e the elementary 
charge, and E(x, t), B(x, t) are the electromagnetic fields which are independent on velocity while being functions of 
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space and time. Finally, V v = d/dv is the velocity gradient operator acting on the phase space distributions. These 



two Vlasov equations ( 18 ) are coupled mutually and to the electromagnetic fields through Maxwell's equations 

(19) 



V x B = /x eo(9E/at) +^oE e ±J dv 3 F ± v, V -B = 

± 

VxE= -(dB/dt), V-E =€o 1 J2e±Jdv 3 F ± 

± 



The second term on the right in the first of these equations is the electric current density; the term on the right in the 
last of these equations is the electric space charge density (divided by the dielectric constant of vacuum, eo). These 
equations already account for the coupling of the field to the particles through the definition of the electric current 
and particle densities as zeroeth and first moments of the one-particle phase space distributions. 

Shocks evolve from infinitesimal wave disturbances; one hence considers two different states of the plasma with the 
physics of both of them contained in the above equations. These two states are first the final average slowly evolving 
state of the fully developed shock, and second the strongly time-dependent evolution of the infinitesimal disturbance 
from the thermal level where it starts up to the formation of the shock. In the first state the shock possesses a 
distinct shock profile while in the second state one deals with initially infinitesimal fluctuations. When the fluctuation 
amplitude approaches the shock strength the two different ways of looking at the shock should ideally lead to the 
same result. According to this distinction one divides all field and plasma quantities, A, into their slowly varying 
averages, (A} 7 and fast fluctuations, <5 A, superimposed on the averages according to the prescription 

A — (A) + 5 A, {5 A) = 

The second part of this prescription breaks down when the fluctuations become very large or non-symmetric, but it 
makes life easier to deal with zero averages of small fluctuations as long as the fluctuation amplitude remains to be 
small. Frequently it is the only way of extracting a solution from the above nonlinear and complex set of equations 
fll8|19D . 



In what follows we will apply different simplifications to all these basic equations referring to the last conditions of 
simplification. Only the last section of the present chapter will, finally, deal with the numerical simulation technique 
which, in fact, will become the most important tool in the investigation of shocks in the remaining three chapters of 
this first part of the book. 



2. Equations for averages and fluctuations 

Let us for simplicity temporarily indicate the averages (• • • ) of the distribution functions and fields by the subscript 
on the unbraced quantities, and the fluctuations by small letters /, e,b. Then, on applying the above prescription 
of averaging to the Vlasov equation, we obtain the kinetic equation for the average distribution functions Fq(v, x, t) 
in the form 

+ v ■ VF± + ^(E + v x B ) • V V F± = _i± <( e + v x b) • V v /±) (20) 
at m± m± 

Here the average quantities are assumed to vary on much longer spatial and temporal scales than the fluctuation 
scales such that the condition of averaging (/, e,b) = remains valid. This average Vlasov equation contains a non- 
vanishing pseudo-collision term on its right which accounts for the effect of the correlations between the fluctuations 



and particles on the average distribution. In contrast to the Vlasov equation, the Maxwell equations (19 1 retain 
their form with the sole difference that the full distribution functions F^ appearing in the expression for the electric 
current density in Ampere's law and in the space charge term in Poisson's equation are to be replaced by their average 
counterparts i 7 ^, yielding 

VxB f j e (aE /^) + /ioEe ± /d U 3 F ± v, V • B = 



V x E^{dB /dt), V ■ E - e^ 1 E e± / dv 3 F± 



(21) 



o 

± 



In order to obtain equations for the fluctuations one subtracts the set of averaged equations from the full set of 
equations and orders the terms for the fluctuation quantities / = F—Fq, e = E— Eo, b = B^Bo- This procedure leaves 
the Maxwell equations unchanged when all quantities appearing in them are replaced by the fluctuating quantities, 
and the fluctuating Vlasov equation becomes 

^r* ■ V/ ± + — (E + v x B ) • V v /± = - — (e + v x b) • V V F± 
at m± m± 

6± -(e + v x b) • V v /± + ^ ((e + v x b) ■ V v /±) (22) 
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Up to this stage the fluctuations are allowed to have arbitrarily large amplitudes; it is only their scales which must 
be much shorter than the scales of the average field quantities. This means for instance that the width of the shock 
transition regions should be much larger than the wavelengths of the fluctuations. 

The last equation is in fact the equation that describes the evolution of fluctuations, however the coupling to the 
average quantities is still so strong in this equation that it can be solved only together with the average equation. In 
particular the average "collision term" appearing on its right provides the greatest complications. It will therefore be 
simplified considerably in treating real problems. 

On the other hand, the "collision term" in the average equation is the term that is responsible for anomalous 
dissipation and is thus the most interesting term in any theory that deals with the evolution of shock waves. For 
a spectrum of properly chosen fluctuations this term prevents large amplitude waves from indefinite steepening and 
breaking and provides the required dissipation of kinetic energy, entropy generation, and shock stabilisation. In its 
general version given above it should contain the whole physics of the shock including the complete collective processes 
which occur before real particle collisions come into play. 

However, the complexity of these equations is still too large for solving them. So one needs further simplifications 
in order to infer about the behaviour of shocks. The simplest and at the same time very effective simplification is 
to ask for the macroscopic conservation laws and the conditions of change of the plasma quantities across the shock 
transition layer which are in accord with the above fundamental kinetic equations. These are the magnetogasdynamic 
equations and the Rankine-Hugoniot jump relations. 



3. Conservation laws 



Following the philosophy of simplification we will first, before asking for the internal processes taking place in the 
shock transition, the generation of dissipation, particle reflection, entropy production etc., look into the global- i.e. 
large-scale - structure of a shock. In order to do this we need consider only the global plasma and field quantities, 
density (N) = N , flow velocity (V) = V respectively momentum density (JW) = A^qVq, pressure (P) = Po, 
magnetic field (B) = B , electric field (E) = E , current density (j) = j , entropy S and so on. 

These quantities are all averages or result from average moments over the global distribution function (F) = Fq. 
Since we will be dealing in the following only with average moments we suppress both the angular brackets and index 
0. The prescription of taking moment of order i is 

M l = / duV.F (23) 



where v l = v. . .v is understood as the z-fold dyadic product. The first three moments are N — J dv 3 F, NV = 
J dv 3 vF, P = m J du 3 (v — V)(v — V)F. Clearly, the diagonal of the pressure tensor P gives the average energy 
density and also defines the local temperatures Tu , Tj_ parallel and perpendicular to the average magnetic field. 



Operating in the usual way with these definitions on the average Vlasov equation ( 20 ) produces the well-known full 
- i.e. infinite - set of magnetogasdynamic equations for the infinite chain of moments of F for each particle species 
± = e,i. The first two of them are 

^-+V-(JW)±=0 (24) 
at 

d (iVV) + V • (NVVf + — V • P ± =^±^(E + V ± x B) + / d« 3 vC ± (25) 
at m± m± J 



where by C the pseudo-collision term on the right of Eq. (20) is meant. Because this term conserves particle number 



(or mass) the zeroeth moment of it vanishes identically and does not contribute to the first (zero-order) of the above 
moment equations. Wave particle interaction neither changes particle number nor mass density. These are strictly 
conserved as is shown by the above particle umber conservation equation. 

In the first order moment equation it produces a wave friction term that has the explicit form 



1 / 1 9 I K\ 

\ — ^r( e x b) 

m± [ /ig at 



e /_2\ , 1 /i_2\A i ( _ /__\ , 1 



+ — (b^) I- eo (ee) + -(bb; 
1 / \ Mo 



(26) 



All these terms arc in fact of nature ponderomotive force-density terms contributed by the average wave pressure 
gradients; the first term results from the wave Poynting moment, the second is the gradient of a pure isotropic wave 
pressure, the third is related to wave pressure anisotropy. The inverse proportionality of this entire expression to the 
mass shows that the main contribution is due to the electron momentum density equation. The effect on the ions 
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FIG. 7: 



A sketch of the planar thin (width A) shock geometry with n shock normal, upstream and downstream regions, bilk flow 

velocities and magnetic field vectors. 



can be neglected as they are (in the non-relativistic case considered here) insensitive to ponderomotive effects. We 
note in passing, that it is this term which while affecting the motion of the electron gas will be responsible for the 
appearance of anomalous collisions, anomalous resistivity and viscosity, which we will discuss at a later occasion. 

The two above equations do not form a complete system of equations. The first contains number density flux, the 
main constituent of the second equation which, as a new entity, contains the pressure. For P one, in principle, can derive 
an energy conservation (heat conduction) equation which would contain the new quantity of heat flux, the next higher 
moment. On the other hand, one can replace the pressure equation that follows from the energy conservation law, 
by equations of state, P(N, 7, Ty , T±), which express the pressure tensor components through density, temperature, 
adiabatic coefficient 7 etc. This is the usual procedure applied when investigating shock solutions. One should, 
however, be aware of the fact that equations of state in non-equilibrium are merely approximations which hold under 
certain conditions of cither isothermality - which does not apply to shocks as they are not in thermal equilibrium - 
or adiabaticity. The latter condition is quite reasonable in dealing with the fast processes taking part in the shock 
environment when the flow passes across the shock front in a time so short that thermalization becomes impossible. 

The idea is to apply the momentum equations to an extended shock that represents a ramp in real space. In the 
spirit of our discussion in Chapter 1, the first step to do this is assuming that the shock is a thin planar discontinuity 
that moves at a certain shock velocity U in the shock normal direction n across the plasma. If we confine all the 
micro-processes to the interior of the shock plane, i.e. if we go far enough away from the shock plane upstream and 
downstream, then we can apply the above dissipation- free average conservation laws to the shock and ask only for the 
differences in the plasma and field parameters between downstream and upstream of the shock, trying to express the 
downstream values in terms of the undisturbed upstream flow and field values. In doing this, we completely neglect 
the "pseudo-collision" terms on the right of these equations, since all physics that is going on will be confined to the 
transition region as wide as it can be. For a plane rigid stationary shock surface this assumption is good enough. 
However, when doing so, with the above separate conservation laws for electrons and ions, we immediately run into 
severe problems even in the simplest completely interaction-free case. The reason is that electrons and ions because 
of their different mass behave completely differently while at the same time cannot be treated separately as they are 
coupled through charge conservation and electrical neutrality and through their unequal contributions to the electric 
current density and therefore to the fields, a difficulty that has been discussed by [57] . 
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III. RANKINE-HUGONIOT RELATIONS 



In order to overcome this difficulty one is forced to further simplification of the conservation equations by adding up 
the electron and ion equations [cf., e.g. , [2J as for one of the many accounts available in the literature]. To this purpose 
one must define new centre-of-mass variables 



± 



1 



TV = 



E± m ± 



v = 



E± m ± N± 



This leads to the magnetogasdynamic or MHD equations for a single-fluid plasma. Since the equation of continuity 
remains unchanged from Eq. ( 24 1 it suffices to write down the momentum conservation equation 



d(mNV) 
dt 



V ■ (mTVVV) = —V ■ P + pE + j x B 



(28) 



where P = P e + is the total pressure tensor, and p is the electric charge density p = e(Ni — N e ) which in quasi- 
neutral plasmas outside the shock is assumed to be zero such that the second term on the right containing the average 
electric field vanishes outside the shock ramp. The last term in this equation is the Lorentz force written in terms of 
the average current from Ampere's law 



Moj = V x B 



(29) 



The displacement current can be safely neglected because these equations hold only for very slow variations with 
frequency smaller than the ion cyclotron frequency lu <C u C i — eB/rrii, scales much larger than the ion gyro-radius 
L r c i = Vi±/uj c i, and wave speeds much less than the speed of light. Note that this equation is completely 
collisionlcss. If we would have retained the pseudo-collision term on the right in the electron equation this would 
simply have added an electron ponderomotive force term on the right. 

In fact, together with Maxwell's equations these equations are not yet complete in a double sense. They have to 
be completed with appropriate equations of state for the pressure components, as has been mentioned above, and 
they have to be completed by a relation between the current and the average electric field that appears in Maxwell's 
equations, i.e. with an appropriate Ohm's law. This is found by subtracting the electron and ion momentum 
conservation equations and turns out to be quite complicated [cf., e.g., 3 1 J . In slightly simplified form Ohm's law 
reads 



E + VxB=4fjxB 

eN 



1 



V-P P 



(30) 



eN ' ' " ' e 2 N dt 

Note that on the right only electron terms are contained in this expression. Also, an electron ponderomotive term - 
responsible for anomalous transport effects - would appear on the right if we would retain the pseudo-collision term. 

However, even in this form even though the system is non-collisional Ohm's law is still too complex for treating 
the conservation laws at a shock transition. The reason is that the right-hand side introduces second order spatial 
derivatives into Faraday's law through the pressure gradient and current expressions. One therefore argues that for 
sufficiently flat shock transitions the terms on the right can be neglected. This argument implies that one must go far 
enough away from the shock into a region where shock excited turbulence has decayed away in order to apply global 
conservation laws to the shock transition, this can be done when only the left-hand side in Ohm's law is retained and 
the ideal MHD frozen-in condition holds: 

E=-VxB (31) 

Assuming that the shock is plane and narrow as shown in Figure[7 such that any variations along the shock can 
be ignored and the sole variation is along the shock normal, Eqs. ( 28pT l, the continuity equation and Maxwell's 
equations become all one-dimensional and can be integrated along n over the shock transition (with regular boundary 
conditions at x = ±oo). Applying the definition of the shock normal (??) in Chapter 1 and the prescription for the 
V-operator in Eq. (??) transforms these equations into a nonlinear algebraic system of equations for the jumps [. . . ] 
of the field quantities 



n • [miVVV] + n 



P 



2p 



n ■ [NV] 

— n [BB] 

Mo 

[n x V x B] 
n x [B] 



mNn ■ V 
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mN 
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= 

= 

= 

= 

= 



(32) 
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Here, for simplicity, the pressure has been assumed isotropic. The last equation is the rewritten energy conservation 
equation, where w = c v P/kBN is the ideal gas enthalpy density, c v the specific heat, and ks Boltzmann's constant. 
This system of equations is the implicit form of the Rankine-Hugoniot conservation equations in ideal magnetogasdy- 
namics (ideal MHD). In this version it contains all ideal MHD discontinuities of which shock waves are a subclass, the 
class of solutions with a finite flow across the discontinuity, compressions (in density) , and increases in temperature 
T, pressure P, and entropy S across the discontinuity in the transition from upstream to downstream. 



A. Explicit MHD shock solutions 



We are not interested in the full set of solutions of the above system of jump conditions (32 1. We rather look for 



genuine shock conditions. This requires finite mass flux T = NV n across the shock in the normal direction. The first 



of the Rankine-Hugoniot relations (32) tells that the jump [T] = 0. Hence T = const, and we must sort for solutions 
with T ^ 0, or NiV n i — N 2 V n2 , m order to be dealing with a shock. 

Introducing the specific volume V = (m-ZV)^ 1 , the whole system of jump conditions can be factorized [5] and can 
be written in the form 

In the one-fluid approximation magnetogasdynamic shock waves with T ^ are contained in the expression in curly 
braces which still depends on the jumps in pressure [P] and specific volume [V] and thus on the energy conservation 
equation respectively the equation of state. We will not discuss this equation further as in the following more insight 
can be gained from explicit consideration of a few particular cases. 

Under the special condition that the flow in the upstream Region 1 is along x (anti-parallel to n) and the upstream 



magnetic field Hi = (B± x , 0, B>i z ) is in the xz-plane, and assuming T ^ 0, the jump conditions Eqs. ( 32 1 simplify. 
Since V\,B\ X ,B\ Z ,P\ are known quantities, it is convenient to introduce normalised variables for the corresponding 
downstream values 



^^No ^^V 2 Tl ' 2 ->T 12 ^ 

Ni ' Vi ' mV 2 /2 1,2 ' y/fiomNtf 



where the temperature is taken in energy units. Instead of it we may also use the corresponding thermal speeds 
which by the above normalisation are normalised to V\ . This then yields the following normalised Rankine-Hugoniot 
relations, in which B n = const as a consequence of the vanishing divergence of the magnetic field, 

N 2 V n2 = 1 

V n 2B z 2 — V z 2B n = B z i 

B Z 2B n — V Z 2 = B z \B n (34) 
2N 2 (vi + V t 2 l2 )+B 2 z2 = 2(l + t,?) + J B» 1 
Vn2 + Vz2 + 2B z2 B zl + hv\ - l + 2B 2 zl + bv\ 

The energy conservation equation yields the last in these expressions. There the enthalpy is taken into account giving 
a factor 5 in front of the thermal velocities. These five equations can be combined into a third-order equation for one 
of the downstream unknown quantities, for instance V n 2, expressed in terms of the upstream values 

a 3 V* 2 + a 2 V 2 2 + ai V n2 + a = (35) 

where a = -B 2 n [b 2 zl + B*(l + 5v 2 )}, oi = 25^(1 + 2B? + 5v 2 ) - \B\ X , -a 2 = l + 5vf + SB 2 + %B 2 1: a 3 = 4. Below 
we discuss a few simple illustrative solutions of this equation. 



B. Perpendicular shocks 



For strictly perpendicular shocks we have B n = 0,B z i — Bi,a = 0,ai = —\B 2 ,—a 2 = 1 + bv\ + ^B zl ,a 3 = 4. 
Equation ( 35 ) turns into a quadratic equation yielding the solution 

Vn2 = l\ 1 + ( 1 + ^ B 2 + ! I 



1 + 2 & ) B f 



2Bi 



1 

N2 



1 

^2 



(36) 
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Since the condition for a shock to exist is that the normal velocity V n 2 < 1 in Region 2 we immediately conclude 
that in a perpendicular shock the density and tangential magnetic field components in Region 2 increase by the same 
fraction as the normal velocity drops, and this fraction is determined by the plasma-/? ratio 0i = 2hqNi T\JB\ in 
Region 1, where T\ = T e \ + Tn is the total temperature. The condition on V n 2 implies that the Mach number takes 
the form (now in physical units) 



1< M = 



Vi 



l + 5/?i/6Vki 1 + 5A/6 



(37) 



Here Ma is the Alfven-Mach number which is the flow to Alfven velocity ratio. In cold plasmas or plasmas containing 
strong magnetic fields j3 <C 1, and the Mach number is simply the Alfven-Mach number. Conversely, in hot plasmas 
the Mach number becomes about the ordinary gasdynamic Mach number. 
For the increase in normalised temperature one finds accordingly 



^ = 1 + ^ fi [l-N^ + 2M- A 2 {l-N 2 )} >1 



(38) 



This is always larger than one. Perpendicular shocks cause plasma heating during shock transition time and thus 
cause also increase in entropy 



AS oc In 



1Y\ -v- 1 j_ 
tJ n 2 



(39) 



which holds under the ideal gas assumption. 



C. Parallel shocks 



This case is not well treated in magnetogasdynamics conditions as we have explained earlier. Since the magnetic field 
is normal to the shock it is theoretically unaffected by the presence of the shock which therefore should become purely 
gasdynamic. In the above perpendicular shock jump conditions one can for this case simply delete the magnetic 
terms. However, this does not cover the real physics involved into parallel shocks which must be treated on the basis 
of kinetic theory and with the simulation tool at hand. 



D. High Mach numbers 

This limit applies when the ram pressure of the flow is very high and exceeds the thermal pressure. Then all terms 
including t> li2 can be neglected. Moreover, one usually also neglects the magnetic field in this case, and the shocks 
become then purely flow determined with V n 2 = N^ 1 ~ B^ 1 — \, suggesting that both, the magnetic field and 
density should not increase by more than a factor of 4. 

In fact, the observations in interplanetary space indeed confirm that all shock that have been observed there are 
weaken than this. However, again, this reasoning does not really apply in plasmas because at very high Mach numbers 
other effects come into play which are connected with the kinetic nature of a plasma, electrodynamic effects, and the 
differences in electron and ion motion. Ultimately relativistic effects must be taken into account. These become 
susceptible first for electrons, increasing their mass but at the same time distinguishing them even stronger from the 
inert ions, because the electron dynamics changes completely in the relativistic domain. 

In addition, high Mach number shocks are supercritical and even though it seems that one could treat them in the 
simple way as has been done here, the kinetic effects involved into their physics inhibit such a simplistic interpretation 
of high-Mach number shocks. High-Mach number shocks readily become turbulent, exciting various kinds of waves 
which grow to large amplitudes and completely modify the environment of the shock which cannot be treated any more 
as quiet. In such a turbulent environment shocks assume intermittent character losing stationarity or even identity 
as a single ramp which the flow has to surpass when going from Region 1 to Region 2. Occasionally the distinction 
between two regions only may become obsolete. There might be more than one transition regions, subshocks form, 
the ramp will evolve its own structures. And these structures come and go, are temporarily created and damp away 
to make space for the evolution of other new structures. Probably, high Mach number shocks exist only temporarily 
at one and the same spatial location. They are highly dynamical, changing their nature, structure, shape, steepness 
and intensity along the surface of the shock such that they strongly deviate from one-dimensionality and even from 
two-dimensionality. They arc time-dependent, reforming themselves continuously in different regions of space and 
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thus cannot be described by a simple plane shock geometry of the kind we have assumed. Later in this book at the 
appropriate place we will consider moderately high Mach number shocks when dealing with the extended class of 
supercritical shocks. 

1. Oblique shocks 

Real shocks do not belong neither to the very particular classes of parallel nor perpendicular shocks. Real shocks 
are oblique in the sense that the upstream magnetic fields Bi are inclined with respect to the shock normal n. As 
mentioned earlier, one distinguishes between quasi-parallel and quasi-perpendicular shocks depending on the shock 
normal angle @Bn being closer to 0° or 90°. Since we will treat the properties of these shocks separately in some 
following chapters, we are not going to discuss them at this place. 

Shocks around finite size obstacles will never be really plane. The best approximations to plane shocks are inter- 
planetary shocks. Bow shocks in front of magnetised planets, comets or other bodies are always curved. They assume 
all kinds of shock properties along their surfaces reaching from perpendicular to oblique and parallel. Curving their 
surface implies that the shock normal changes its angle with respect to the direction of the upstream flow Vi. 

Since shocks evolve for Mach numbers A4 ms — V\/c ms > 1 they occupy a finite volume of space only. For 
shock formation the Mach number based on V\ n normal to the shock is relevant. Defining the flow normal angle 
cos6y„ = |Vi-n|, formation of a bow shock in front of a finite size object like a magnetosphere is restricted to flow 
normal angles 

e Vn <cos- 1 (A^- 1 J (40) 
For instance, at a nominal Mach number M ms ~ 8 shocks exist for angles 6y„ < 82.8°. 

IV. WAVES AND INSTABILITIES 

It has been mentioned several times that shocks evolve from waves mainly through nonlinear wave steepening and the 
onset of dissipation and dispersion. Moreover, it is the various modes of waves that are responsible for the generation 
of anomalous dissipation, shock ramp broadening, generation of turbulence in the shock environment and shock ramp 
itself, as well as for particle acceleration, shock particle reflection and the successive effects. The idea is that in a 
plasma that consists of electrodynamically active particles the excitation of the various plasma wave modes in the 
electromagnetic field as collective effects is the easiest way of energy distribution and transport. There is very little 
momentum needed in order to accelerate a wave, even though many particles are involved in the excitation and 
propagation of the wave, much less momentum than accelerating a substantial number of particles to medium energy. 
Therefore any more profound understanding of shock processes cannot avoid bothering with waves, instabilities, wave 
excitation and wave particle interaction. 

A. Dispersion relation 

Waves are a very general phenomenon of most media. However, they do not fall from sky. Instead, they evolve 
from small thermal fluctuations in the medium. Such fluctuations are unavoidable. In order for a wave to propagate 
in the medium a number of conditions need to be satisfied, however. The first is that the medium allows for a 
particular range of frequencies u) and wave-vectors k to exist in the medium; i.e. it allows for eigen-oscillations or 
cigen-modes. These ranges are specified by the dispersion relation D(oj, k, . . . ) — which formulates the condition that 
the dynamical equations of the medium possess small-amplitude solutions. This dispersion relation is usually derived 
in the linear infinitcsimally small amplitude approximation. However, nonlinear dispersion relations can sometimes 
also be formulated in which case D depends on the amplitude as well. 

Plasmas arc elcctromagnetically active media with the electromagnetic field governed by Maxwell's equations. Since 
there the plasma properties enter only through the material equations (i.e. current density j, space charge p), the 
dispersion relation is most easily obtained from them. Moreover both, the current and space charge in a plasma, 
depend on the number densities in the plasma; i.e. the space charge variation can be included into the current 
variation. It is then simple matter to derive the general electromagnetic wave equation for the fluctuating fields, e, 
on a much slower evolving background, (B, E), 

V 2 e - V( V • e) - 60^07^ = Mo ^ (41) 
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The magnetic fluctuation field b is completely determined from Maxwell's equations, and the current appearing on 
the right is expressed conveniently through the space-time dependence of the fluctuation-conductivity tensor <r(x, t) 
as 



j((x,i) 



J dx' J dt'a(x-x',t-t') 



(42) 



an expression that implicitly accounts for causality due to the integration over the entire past of the current up to the 
observation time t. Since the fluctuation current is a functional of the complete set of particle distribution functions 
(through the zeroeth, N, and first, NV, moments) the complete evolution of the fluctuations up to a large amplitude 
shock is contained in these expressions. However, for practical purposes one linearizes this equation by assuming 
that the fluctuation-conductivity tensor is, to first order, independent of the fields e,b. The above equation (41 1 



becomes linear under this assumption and can be Fourier analysed, with wave vector k and frequency u>, yielding 
(with c 2 = l//ioeo the square of the speed of light) 



k 



I — kk — iujfio<r(uj, k) 



•e(w,k) = 



(43) 



The quantities in this expression satisfy the following symmetry relations e(— lu, — k) = e*(uj, k), a{— to, — k) = cr*(u>, k). 
Setting the expression in brackets to zero yields the equation for the linear eigenmodes u> — w(k). For convenience 
we define the dielectric tensor 



e(w,k) = l 



UJ£ 



-a(w,k) 



(44) 



which satisfies the same conditions as the fluctuation conductivity. Then we can finally write the general dispersion 
relation in the compact form as the determinant of the bracketed expression 



D(oj,k) = Det 



I -kk + e(w,k) 







(45) 



The particular linear physics of the plasma is contained in the dielectric tensor through the conductivity. (We note 
in passing that for any classical medium the above conditions together with the dispersion relation are the equivalent 
of the well-known Kramers-Kronig relations of causal fluctuations in quantum mechanics [see, e.g., [32].) In addition 
to the solution of the dispersion equation ( 45 1 the first problem consists in the determination of the linear fluctuation 



conductivity tensor which enters the dielectric tensor (44 1. For this one needs to go to the appropriate plasma model. 



However, we repeat that without this linear step one cannot obtain any susceptible information about the nature of 
a shock wave. This we have explained in breadth in Chapter 1 and the preceding sections. 

The linear dispersion relation Eq. (45 1 has plane wave mode solutions of the form cx expi(k • x — cut), which are 



eigenmodes of the particular plasma model which is described by the kinetic equations (or appropriate simplifications 
of the kinetic equations) of the plasma. The dynamics of the plasma enters through the wave conductivity tensor 
which can be determined from the Fourier transformed expression for the current density 



j(k l w)=J3±e± J d« 3 v/ ± (v ) k,w) = <r(k,a;)-e(k J w) 



(46) 



as the first moment of the fluctuating part / of the distribution function F ± . Solving for the integral and expressing 
f ± through the electric wave field e yields th e wa nted form of the wave conductivity tensor a. The problem is thus 



reduced to the determination of / from Eq. ( 22 1 where we drop the average terms on the right-hand side retaining 
only terms linear in the fluctuations: 



v-V/ d 



(E + v x B ) • V v / d 



(e + v x b) • Vv^ 



dt m± m± 

Operating with a Fourier transform on this equation then yields the following expression 

,± _ ie± i(e + v x b) ■ V V F^ 



1 



ie± 



m±(k • v — u>) 



(E + v x B ) ■ V v 



P 



m± 



(47) 



(48) 



which determines f ± in terms of the average and fluctuating field quantities, which is just what we want. One can 
now make assumptions about the average fields and distribution function in order to explicitly calculate f ± . Usually 
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FIG. 8: Left: Two different types of dispersions in the (real) {u> r , k)-plane. Short waves with concave dispersion have slower group 
velocity than long waves and form a trail on the long wave. Short waves with convex dispersion move faster than long waves. Right: The 
effect of this difference in sidewave velocity on a laminar subcritical shock wave. Convex dispersions produce run-away waves which 
appear as spatially damped upstream oscillations (in the B-field, for instance). The trailing waves resulting from concave dispersion 
occur as downstream spatially damped oscillations. Maximum wave amplitudes are observed near the shock ramp in both cases. 

these assumptions are Eo = 0,Bo = B^i with Bq = const. Then the operator V v = —zd/d<fi on the left simplifies 
to a mere derivative with respect to the gyration angle (f> of the particles. Further assumptions on F are that the 
average distributions are gyro-tropic, in which case the integration with respect to <fi becomes trivial. 

With such assumptions it is not difficult so still tedious to solve for and finally get the conductivity tensor. The 
respective expressions have been given in various places and will not repeated here. Good references among others 
are [ID] , [ST] , [5] . One can even include weak inhomogeneities in the average distribution function and fields which is 
needed when considering an inhomogcncous initial state like a given soliton or shock structure and investigating its 
prospective stability. In this case the plasma background is not homogeneous anymore because a soliton or shock has 
already evolved in it and has locally modified the plasma. Any wave modes which will be excited on this modified 
background will then not be influenced only by sideband formation, steepening, nonlincarity and dispersion but also 
by the change of the plasma properties from location to location. This implies that the waves themselves change 
character and properties across a shock. 



1. Damping/growth rate 



The solutions of the dispersion relation arc in most cases complex, and for real wave vector k can be written as 
c<j(k) = ca r (k) + i7(uv,k), where the index r indicates the real part, and 7 is the imaginary part of the frequency 
which itself is a function of the real frequency and wave number, because each mode of given frequency can behave 
differently in time, and the wave under normal conditions will be dispersive, i.e. it will not be a linear function of wave 
number. In most cases the amplitude of a given wave will change slowly in time, which means that the imaginary 
part of the frequency is small compared to the real frequency If this is granted, then 7 can be determined by a simple 
procedure directly from the dispersion relation Z?(w,k) = _D r (a;, k) + iDj(o;,k), which can be written as the sum of 
its real D r and imaginary Di parts because a small imaginary part 7 in the frequency changes the dispersion relation 
only weakly, and it can be expanded with respect to this imaginary part. Up to first order in 7/0; one then obtains 

A-K,k) = 0, 7 K,k) = - (49) 

The first of these expressions determines the real frequency as function of wave number uj r (k) which can be calculated 
directly from the real part of the dispersion relation. The second equation is a prescription to determine the imaginary 
part of the frequency, i.e. the damping or growth rate of the wave. 
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2. Remarks 



Two remarks on the dispersion relation are in place. First, the weak damping/growth rate solutions u) — w(k) of 
the above general dispersion relation - themselves called dispersion relations - arc also of use in the weakly nonlinear 
case. They can be understood as the lowest-order expansion term of a more general nonlinear dispersion relation 
lu = w(k, |e,b| 2 ) which depends weakly on the wave amplitude or wave energy |e,b| 2 . When taking into account 
higher order expansion terms in the wave amplitude it produces other non-linear equations which govern the amplitude 
evolution of the wave under consideration. Such an equation is the non-linear Schrodinger equation we will get familiar 
with when discussing transport processes. 

Second, from the dispersion relation uj r — Lu r (k) one can infer in which way steepening of a wave is compensated 
by the dispersion of the wave. Figure 8] on its left shows two typical cases of (real) dispersion curves of low frequency 
waves in the (u> r , /c)-plane from which shock waves could evolve [after Ho] . Both curves have in common that they 
exhibit linear dispersion at long wavelengths, i.e. at small wave-numbers k, with slope giving the phase velocities of 
the waves. In this region all nonlinearly generated sidebands have same phase and group velocities causing broadening 
of the wave spectrum and steepening. However, at higher wave-numbers the dispersion curves start diverging from 
linear slope, one of the waves turning convex, the other concave. These turnovers imply a change in phase and group 
velocities. The convex dispersion implies that shorter wavelengths generated in the convex part of the dispersion 
curve move faster than the long waves. They will thus catch up with the long wavelength wave and run away ahead 
of the wave forming upstream precursors of the wave as shown for the shock in the lower part on the right. On 
the other hand, for the concave dispersion shorter wavelength waves fall behind the long waves. They represent a 
wave trail following the large amplitude long wave as is shown for the shock in the upper part on the right. Hence a 
simple glance at the dispersion curves already unmantles the possible properties of the expected nonlinearity and the 
structure of the shock. 

A word of caution is in place here, however. This reasoning does not hold for all shocks but for subcritical laminar 
shocks only. Supercritical higher Mach number shocks will behave in a more complicated way being much less 
dependent on dissipation and dispersion. 



The waves from which a shock forms are the lowest-frequency plasma modes that are excited under the particular 
conditions of the shocked plasma. In low-/? (cold) plasma these are the tree fundamental MHD modes. Since all 
collisionless shocks in the heliosphere are magnetised the magnetic field has to be included, and the global shocks are 
not purely electrostatic even though subshocks developing in them can well behave approximately electrostatic. 

We arc already familiar with the three low-/? magnetogasdynamic modes, the fast and slow magnetosonic and the 
intermediate Alfven waves. These are the lowest frequency eigenmodes of a homogeneous not necessarily isotropic 
plasma, i.e. when a small disturbance is present in the plasma it will propagate in one or all of these modes. Their 
dispersion relation follows from 



which depends on the parallel and perpendicular components of k only in a very simple way. Moreover, it is a purely 
real dispersion relation lacking any imaginary part and therefore also any damping which is of course typical for a 
low frequency dissipation-free plasma. It is a different question of how these modes can be excited, and we will come 
to this at a later stage. 

The phase velocity c ms of these modes has been dealt with already in Eq. (??). Since these waves are linear waves 
with no dispersion, their dispersion relation is simply u> — kc ms {6) with 6 the angle between the wave vector k and 
magnetic field B. Figure|9] shows the real space angular dependence of these three phase velocities for two special 
cases. Clearly in the direction perpendicular to the magnetic field only the fast mode propagates and, hence, strictly 
perpendicular MHD shocks are fast shocks as has been noted. In the direction parallel to B all three waves can 
propagate. 



B. The MHD modes - low-/? shocks 




(50) 
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FIG. 9: Wave vector diagram of two cases of MHD waves in the plane of the magnetic field. 



1. Magnetosonic solitons in cold plasma 



We will now show that the method of the Sagdeev pseudo-potential can be used to understand the formation of a fast 
mode solitary wave (or soliton) propagating strictly perpendicular to the magnetic field B. This has first been shown 
by [9]. What results from this procedure will not yet be a shock, because, as we have noted before, shock formation 
requires the presence of some kind of dissipation, while the equations on which the present theory is based are strictly 
dissipation-free (remember that we have dropped the correlation terms on the right-hand side of the kinetic equation 
before deriving the moment equations and that we have not yet discussed any way of how dissipation occurs when 
these terms are taken into account). 

In order to find the stationary solutions we are looking for, one must retain the nonlinearity in the stationary one- 
dimensional quasi-neutral magnetogasdynamic equations. This nonlinearity appears in the convective term V n V n V n 
in the equation of motion. From constancy of normal flux [T] = one has NV n = NiV\ and E y — N\Bi where the 
index 1 means undisturbed values far upstream. Since only electrons contribute to the current by their drift in the 
crossed electric and magnetic fields E = —V xB,B = Bi in the shock frame, we must retain a small component 
E n = —BV y across the shock, effectively produced by the difference in electron and ion motion and causing the shock 
current to flow in y-direction in the shock transition. It occurs in the stationary equation of motion on the scale of 
the shock transition, i.e. on the scale of the ion gyro-radius, and as current j v = —eNV y y in Ampere's law 



mNV n ^- = -eNBV y , ^ = -fi eNV y (51) 



Combining these equations yields the normal fluid velocity as function of the magnetic field and the initial bulk flow 
velocity V\ at infinity 



2 fM) mN 1 V?) X \ 2 V 2 



Vn = V 1 1- ~ r* )=V X (52) 



showing that the bulk velocity decreases from V\ when the magnetic field B increases. The second term in the first 
parentheses is the difference in the ratio of magnetic pressures B 2 /2fiQ at the location under observance and B 2 /2^ 
at infinity upstream to the kinetic pressure mNV 2 at infinity. This is written in terms of Alfven velocities in the 
second parentheses on the right. 

With the help of this expression, E y and Ampere's law we obtain 



B 2 - B\ 
2fi mN 1 Vi 



- 1 



f < 53 > 

da; 
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which can be used to eliminate the velocity components and obtain an equation for the variation in the magnetic field 
given in the form of the energy conservation equation of a pseudo-particle of mass 1 and velocity B in the Sagdeev 
pseudo-potential S(B): 



Hd^J = ^ ( ) — 2Xl[l-(B^-Bj)/2^mN 1 V 1 ^ (54) 



where A e = c/uj pe is the electron inertial length (electron skin depth). The electron inertial length is the only length 
scale that appears in the above equation when we perform a dimensional analysis. It therefore turns out that the 
characteristic width, A ~ \ e , of the magnetosonic solitons is of the order of the electron skin depth. This can 
also be seen when for small amplitude disturbances the last expression is expanded with respect to B. Defining 
b= (B - B 1 )/B 1 ,b m = (B m - Bx)/B u £= x/A yields to first order 

— ~±b(b m -b)' -» b -TTT fti^' A =^T~ 55 

d£ (1+cxp-ICI) 2 sfb m 

giving the above scaling of the soliton width. In addition the inverse scaling of A with the maximum soliton amplitude 
b m is reproduced. 

As discussed before, S < is required for solutions to exist. It is clear that in the absence of dissipation no shock can 
emerge from these stationary waves. They are solitary waves, stationary wave structures of finite spatial extensions 
and amplitude. 

Returning to the original variables, the maximum soliton amplitude is obtained from S(B m ) = as B m = 

yields that B m < 3Bi, and consequently the Alfvenic Mach num- 
oser the Mach number approaches the maximum Mach number 2 the 



B 1 (2V 1 /V A i - 1) which together with Eq. ((52 
ber for the solitons to exist is M. a < 2. The c 



narrower the solitons become. This means that they steepen and, for Mach numbers M. > 2, will overturn and break, 
because the dispersion does no longer balance the nonlincarity. Solitons cannot exist anymore at Mach numbers such 
high. 

When a way can be found to generate dissipation in the region occupied by the soliton, then a magnetosonic soliton 
can evolve into a shock wave. Sagdeev's idea was that this can happen when the soliton becomes large amplitude and 
narrow enough such that in the steep rise of its crest sufficient dissipation could be generated by anomalous collisions 
and anomalous friction. These anomalous collisions would generate sufficient entropy that the states on the two sides 
of the soliton would differ from each other and the flow across the soliton would irreversibly change. In this case the 
soliton would turn into a dissipative subcritical laminar shock. Such shocks will be discussed in Chapter 3. Formally 
we may, of course, model the dissipative effect by simply defining some "collision frequency v" and introducing a 



"collision term" on the right of Eq. ( 54 ) . In order to do this we go one step back to the equation from which ( 54 1 has 



been obtained and add the collision term there: 

d 2 B dS(B) dB 

d^ = -^5--^ (56) 

This equation is modelled exactly after the equation of a damped oscillator, where dB/dx is the velocity. We should 
note here that this modelling has not been justified yet and in fact is not justified by any of our arguments yet. It 
not only requires the proof of the existence of an anomalous collision frequency v, it also requires the proof that from 
the kinetic equations containing the correlation terms an equation of the above structure can be derived. 

Ignoring these objections and exploiting the analogy with the damped oscillator we may conclude from the equation 



( 56 1 for the damped oscillator that the inclusion of anomalous "collisions" will dissipate the kinetic energy of motion 
of the pseudo-particle during its oscillation in the Sagdeev pseudo-potential S(B) until the particle will finally come to 
rest at the bottom of the potential well. This is the case we have discussed earlier in connection with the Korteweg-de 
Vries-Burgers equation. It is drawn schematically in Figure[6] when a shock wave forms from the soliton. The value 
of B at minimum in the Sagdeev pseudo-potential where the pseudo-particle ultimately settles is the magnetic field 
level 

[(8/3! + 1)5 -lj (57) 



Bo — — Bi 
2 



far downstream of the shock that has formed in this dissipation process from the magnetosonic soliton. This value is 
determined by the upstream plasma-/? value. 

Taking this for granted, we can conclude that the damped oscillations the pseudo-particle performs on its damped 
downward path in the Sagdeev pseudo-potential are the spatially damped oscillations of the field B(x) downstream of 
the shock. Moreover, the shock possesses an overshoot in B at shock position B^ < B ov < B m which is smaller than 
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FIG. 10: Left: Dispersion relation of perpendicular magnetosonic waves with concave dispersion at large k, where the wave effectively 
becomes a lower-hybrid wave with frequency close to the lower-hybrid frequency loi^. Right: Dispersion relation in cold plasma at nearly 
parallel propagation. The ion branch starts from the left-handed Alfven modes and goes in resonance at the ion-cyclotron frequency U) c i 
with concave dispersion. The electron branch starts from the right-handed Alfven mode, has convex dispersion, passes through the lower 
hybrid mode into the whistler mode uiih <u« ui ce , assumes concave dispersion, and finally goes in resonance at the electron-cyclotron 

frequency u> C e, where it becomes the electron-cyclotron wave. 



B m but larger than B2 . The existence of damped downstream oscillations is in agreement with the concave shape of 
the dispersion relation of magnetosonic waves which, for large k, become dispersive and approach the lower-hybrid 
branch. This can be seen directly from the dispersion relation for perpendicular (6 = 90°) propagating magnetosonic 
waves 

u ms = V A k±(l + k 2 \ 2 e y i (58) 

For k 2 X 2 e < 1 the wave has constant phase velocity and is non-dispersive, becoming gradually dispersive with increasing 
k when the effective phase velocity decreases. Hence, dujjdk < and the dispersion is concave, for very large k\ e 3> 1 
approaching the lower hybrid frequency ioih = ^/uj ce uj C i — Lu ce y/m e /mi where it flattens out, as shown in Figure 



10 



At oblique angles 90° — <5 ^> \J m e /rrii « 1/43 the dispersion is inverted, and dujjdk > 0. Here the shorter waves run 
the soliton out and appear on the upstream side as spatially damped oscillations. Now their scale is the ion inertial 
length A, = c/uj pi . However, shocks with convex dispersion where the shorter waves outrun the soliton will not exhibit 
a sharp shock profile. Rather they will be oscillating shocks with smoothed out ramp. 

The theory presented above applies to a low-frequency plasma of velocity V± » v e ,Vi larger than the thermal 
velocities of the plasma components. When the relation between V\ and the thermal velocities changes, one must take 
into account thermal effects. These change the nature of the solitons and shock substantially. These changes will be 
discussed in greater depth in Chapter 3. 



C. Whistlers and Alfven shocks 



Many waves are capable of dispersively evolving into solitons or other similar stationary wave pulses if only their 
dispersion relation allows it. Among those waves we here consider only two particular cases, whistlers in cold plasma 
and Alfven waves in low but finite (3 conditions. 
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1. Whistler solitons in cold plasma 



We now return to the cold plasma dispersion relation including both kinds of particles, electrons and ions. For parallel 
propagation k = fc||Z, where B||z, and the angle 9 = 0, the dispersion relations of linear plasma waves are 

k 2 c 2 ui 2 , s 

—a" = 1 - , \ ~ 1 1 \ 59 

0J Z LU(U! =F U) ce ) UI(U!±U! c i) 

The frequencies in the nominators are the electron and ion plasma frequencies, u) pe , P i, respectively. In the denominators 
appear the electron and ion cyclotron frequencies, ui ceci respectively. These relations describe right-hand and left- 
hand polarised waves according to the upper and lower signs. Basically two branches coming from the resonances in 
the denominators are described by this relation. Figurc[lO]on the right shows a plot of the two dispersion curves. For 
the evolution of shocks the most interesting part is the shaded whistler-mode dispersion relation. 

The whistler dispersion relation is the (upper sign) electron part of the above dispersion relation. Neglecting the 
non-resonant ion contribution it reads 

k 2 c 2 uj 2 

— = 1 + - pe - 60 

with the second term being large because of the resonance in the denominator. Thus one can also neglect the 1 on the 
right finding that solutions exist only for u> < u> ce as is also seen in the above drawing. Then the dispersion relation 
becomes 

k 2 \ 2 e ~ Uj/(uj ce - Uj) -» Id ~ UJ ce (l + p^2^) ( 61 ) 

which exhibits its concave character confirming that short wavelengths whistlers will fall behind the main shock pulse. 
Nonlinear analysis of these waves goes back to [33], @B], and [27] and is based on the fluid equations we used before. 
Let the plasma again be moving in x-direction antiparallel to the shock normal n and write the magnetic field in polar 
coordinates as B = B±(0, cos 9, sin 9); then one again obtains the canonical Sagdeev form of the first integral of the 
equation of motion of a pseudo-particle at pseudo-position B± and pseudo- velocity dB±/dx as 

1 fdB ± \ 2 



2 UfJ = ~ s{Bl - ) (62) 



The Sagdeev pseudo-potential is a complicated expression which simplifies considerably for a uniform upstream state. 
We introduce the normalised variables = B±/B± m ,£ — x/X e and [3j_ = B\j fi mNiV 2 writing 

2 [~dj) - _S(6±) -"8^ m (l-Wx) 2 (63) 

Solitons exist for B± < B± m which is the maximum whistler soliton amplitude, and for ft± > |% For the maximum 
amplitude we have 

B ±m < ^nom^V? - ^m,/m e B 1 « 43 B x (64) 
and the whistler soliton velocity V\ ^S> Va\ yielding a soliton Mach number range of 

\/rrii/Am e < Ma < y / m i /2m e — > 22 < Ma < 30 (65) 

which identifies the whistler solitons as being high-Mach number solitons, indeed. In case they evolve into shocks, 
these shocks are high-Mach number as well. This might cause other effects which have not been considered so far. 
Hence the present formal theory must be taken with caution in application to real problems of much lower Mach 
numbers. One of the neglected conditions is quasi-neutrality, which demands that eoE/eNAx <C 1. This leads to the 
further restriction on B± and the pulse width A 

Bx/Bi « {2^c 2 /Vl 1 Y ~ 1.5c/V A1 , A ~ Ax ~ Xi(Bi/B±) (66) 

where = c/uo pi is the ion skin depth (ion inertial length). It follows that these whistler pulses should have quite 
large characteristic widths and, moreover, characteristic frequencies u> ~ ^0Jih(B±/Bi), far below the lower hybrid 
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FIG. 11: Phase velocity diagram of the three Alfven wave modes in the (B, k)-plane of the magnetic field and wave vector for the two 
cases when the alfven velocity is larger left or smaller /right than the electron thermal speed vth.e- The ordinary Alfven wave describes a 

circle in this plane. The two other phase velocities describe deformed curves. 



frequency. These properties identify the whistlers as right-handed (rotating clockwise along x) high-frequency Alfven 
wave pulses. These low frequency whistler/high frequency whistler-Alfven waves can indeed been excited by a cold 
shock-reflected ion beam as will be shown below in the section on ion beam instabilities. 

In view of our remarks we do not give rigourous derivations of these approximate formulae. It will turn out later 
when discussing numerical simulations that whistlers do indeed occur at shocks and have been observed early on as 
well in laboratory experiments on collisionless shocks [e.g. .112) that, however, other fluctuations driven by ion reflection 
are of greater importance in structuring supercritical Mach number shocks. Still it is highly probable that whistlers 
are excited in shock waves as the conditions will be in favour of them when sufficient free energy is available in the 
shock front because of several reasons, one of them electron heating in the perpendicular direction [early observations 
suggested their presence near shocks in space, see, H5] . In this case whistlers become very important for producing 
dissipation via a short- wavelength instability, called decay instability, which had been predicted by [19 and for which 
evidence has been found in the above laboratory observations by [T2"] . 



2. Alfven solitons at finite-f3 



Alfven waves are non-dispersive. However, when the plasma temperature increases, dispersion in the direction per- 
pendicular to the magnetic field sets on. From the general dispersion relation in the very low-frequency limit one then 
obtains [e.g., [2] for the frequency of the Alfven wave 



■„2 „2 
v ± r ci 



where the ion gyro-radius is slightly modified with temperature according to r 2 ci — > r^(|+T e /Tj), and V\ = B 2 / fj, mN 
is the square of the Alfven speed. This dispersion relation describes two kinds of Alfven waves depending on k± ~ r" 1 
or k± r~j X~ . These modes become important when the plasma-/? < 1. The phase velocities of the two modes together 



with the ordinary Alfven wave are shown in Figure 1 1 The modes differ in their dispersive properties from ordinary 
Alfven waves as they propagate oblique to the magnetic field, i.e. the wave energy propagates under an angle to the 
magnetic field for k± is independent of kii . 

At finite temperatures 1 > (3 > m e /m.i the wave is called kinetic Alfven wave. Its perpendicular wavelength becomes 
the order of the ion gyro-radius, and the phase speed parallel to the magnetic field increases. At low temperatures 
[3 < m e /mi the wave is called shear or better inertial Alfven wave with perpendicular wavelength comparable to the 
electron skin depth, and the effective parallel phase velocity decreases. In terms of prospective Alfvenic shocks this 
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means that a kinetic Alfven shock in direction parallel to the external field will support oscillations upstream of the 
shock ramp, while an inertial Alfven shock will support downstream oscillations and thus possess a sharp shock ramp. 

The dispersive properties of both kinetic Alfven modes enable the existence of stationary wave pulses. These have 
been discovered first by [26] for the kinetic mode. For the low Alfven frequencies quasi-neutrality is a good assumption. 
Moreover, the electrons have time enough to behave Boltzmann-like. Hence with the parallel electric potential <j)\\ we 
have 

TV e = TV exp( e </>||/T e ) (68) 

In the perpendicular direction we use the electric potential <pj_ . Maxwell's equations, the nonlinear ion momentum 
conservation, and Poisson's equation then reduce to 



dB± 



ViVf(£i-0||) = MoV||^ (69) 

dt Bouj C i \ 1 dt 

The field-aligned current jy is carried by the hot electron component, such that its divergence is given by V||j|| = 
e(dN e /dt). introducing dimensionless variables £ = xuj C i/ci a ,C, — zaj p i/c,T — u> C it,N' = N/Nq and measuring the 
potentials in units of e/T e , then transforming to a comoving coordinate system r\ = + «n£ — t, the whole system 
of equations is reduced to the nonlinear equation 

«i^'^^ = (l-^)(iV'-Kf) (70) 

which is in the form suited for the Sagdeev pseudo-potential method, yielding (dTV'|d?7) 2 = —S(N', kii , k±) with 

2TV' r i 
S(N',k\\,k ± ) = -^ k 2 [(1 ~ N')(N' + kI) + (1 + kI)N' In N'^ < (71) 

for soliton solutions to exist. These solutions give the density as function of the linear coordinate rj. Interestingly, 
there solutions which are dilutions and solutions which are compressions. The condition for existence of soliton 
solutions is independent of the perpendicular wavenumber. hence it is the parallel electric field that is responsible for 
the formation of solitons and balance of the nonlinear steepening. Solitons form only in parallel direction with the 
magnet field being inclined to the soliton which in the perpendicular direction is flat. When such a soliton attains 
dissipation and turns into a shock, it will be a quasi-parallel shock preceded by damped upstream waves that have 
outrun the shock ramp. That the shock will be quasi-parallel can be easily seen from the fact that k± 3> fcii , the 
shock front will be perpendicular to the external field, and therefore b± <C b» as required for a quasi-parallel shock. 



Inspection of Eq. (71 1 shows that the Sagdeev potential vanishes at TV' — 0, TV' = 1, and N' m . Compressive 
(rarefaction) solitons occur at TV' > 1 (TV' < 1). Only compressive solitons are of interest in shock formation. The 
maximum amplitude N' m of compressive solitons follows from setting the bracket to zero. It is approximately given 
by the solution of TV' + k?< p» (1 + In TV' which, for «jj = 1 is TV r '„ rj 3. The minimum of S(N') for compressive 
solitons is found by taking the derivative of the bracket and putting it to zero. Setting TV' = 1 + n' and expanding 
the logarithm one finds the minimum of the Sagdeev pseudo-potential trough at 

TV^ AS = l + n'« — -J- <3 (72) 

I 

This is the prospective maximum amplitude of a kinetic Alfven shock (KAS) evolving in the presence of dissipation 
from a compressive kinetic Alfven soliton. The compression at a KAS ramp is thus limited to a factor < 3. Finally, the 
value of the Sagdeev pseudo-potential at its absolute minimum gives the steepest gradient of the density in the shock 
ramp, (dTV'/d?7) max = |5 , (TV ( ^ AS ) | at the turning point of the density in the ramp. This ratio provides an estimate of 
the width of the KAS-shock pulse 

Akas ~ Kas/I^as)! > |3 - Kjjr 1 (73) 

in normalised units (where in the final estimate at the right we have used the maximum shock amplitude) . It depends 
on the given values of Ky while being independent on kj_, in agreement with the fact that only values parallel to the 
field will affect the shock structure. 
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We will not investigate the case of inertial Alfven wave solitons. In the heliosphere these shocks occur only in the 
auroral zones of strongly magnetised planets and, possibly, also deep in the solar corona only in region of very strong 
magnetic field and comparably small temperature. The related solitons have been discussed elsewhere [3J EI]- Shocks 
forming there should always be rarefactive, containing very little plasma and having highly oscillating wakes on their 
downstream side containing short wavelength inertial Alfven waves which have been left behind the faster shock ramp. 
In fact such very diluted plasmas have indeed been inferred in the auroral zones of planets like Earth and Jupiter. 



3. Remarks on the generation of dissipation 

|46j gave an idea of how dissipation can be generated in a shock ramp realising that shocks must contain a - large 
scale electric potential drop (or its equivalent as, for instance, an equivalent electric field corresponding to the 
shock ramp density gradient) at which low energy ions will be reflected. (Note that electron will instead become 
accelerated by this potential across the shock.) No matter how few those ions are, they will return into the upstream 
medium, where they (as |46j had noted) become accelerated tangentially along the (perpendicular) shock surface by 
the upstream convection electric field until gaining enough kinetic energy to overcome the shock potential barrier, 
passing the shock, and escaping downstream. These ions form a current in front of the shock that carries free energy 
and will ultimately become unstable with respect to the two-stream instability, scattering upstream electrons and in 
this way cause dissipation. The physical mechanism of this process will be discussed later in this chapter in the section 
on anomalous transport. First we need to be informed about the instabilities and waves that arc relevant with respect 
to shock formation. To these mechanisms we will continuously return when discussing the different types of shocks 
and the corresponding numerical simulations. Without them neither the existence nor the structure of shocks in 
collisionless plasmas can be understood. The models presented so far cannot give more than hints in which direction 
one has to pursue. Shock physics is too complicated for analytical theory. 



D. Instabilities 



An instability is the reaction of an active medium like a warm plasma to the presence of free energy. Since the 
available free energy keeps the plasma away from thermodynamic equilibrium, restoration of equilibrium becomes 
necessary. This is most easily done by exciting fluctuations to amplitudes large enough for either causing dissipation 
or transporting the energy away to a location where it can be dissipated by other processes. As a consequence a 
wave will start growing out of the thermal fluctuation background. The selection of the frequency and wavelength 
range usually is ruled by resonance with the number of particles which carry the free energy. However, other ways of 
exciting waves are also possible. The wave that survives is that with the fastest grow, and least damping. [Theoretical 
overviews of instabilities and dissipation related to shocks have been given by [TH1 HH 113 [5U [53] 

Since most instabilities follow this recipe of growing out of the thermal background they start as infinitesimal 
disturbances which can be described by linear dispersion theory. Hence instability theory in plasma can be based on 
the dispersion relation Eq. (45 1 and for weakly growing waves on the expression for the growth rate (49 1. Plane linear 
wave modes possess a phase factor exp i(k - x — u>t). It is clear that instability sets on whenever the imaginary part 
of the frequency becomes positive; 7(w r , k, . . . ) > is the growth rate of the instability of the particular wave mode 
that becomes unstable. 

In the context of collisionless shocks the instabilities of interest can be divided in two classes. The first class contains 
those waves which can grow themselves to become a shock. It is clear that these waves will be of low frequency and 
comparably large scale because otherwise they would not evolve into a large macroscopic shock. We have already 
discussed a few candidates and their nonlinear evolution in the previous sections, among them magnetosonic, Alfvcnic 
and whistler modes. In this section we will investigate a number of waves which form secondarily after an initial 
seed shock ramp has grown in some way out of one of these wave modes, these are ion modes which have now been 
identified to be responsible for structuring, shaping and reforming the shock. In fact real oblique shocks - which are 
the main class of shocks in interplanetary space and probably in all space and astrophysical objects - cannot survive 
without the presence of these ion waves which can therefore be considered of the wave modes that really produce 
shocks in a process of taking and giving between shock and waves. 

The second class are waves that accompany the shock and provide anomalous transport coefficients like anomalous 
collision frequencies, friction coefficients, heat conductivity and viscosity. These waves are also important for the 
shock as they contribute to entropy generation and dissipation. However, they are not primary in the sense that they 
are not shock-forming waves. Among them there is another group that only carries away energy and information 
from the shock. These are high-frequency waves, mostly electrostatic in nature, produced by electrons, or when 
electromagnetic they are in the free-space radiation modes. In the latter case they carry the information from remote 
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objects as radiation in various modes, radio or x-ray to Earth, informing of the existence of a shock. In interplanetary 
space it is only radio waves which fall into this group as the radiation measure of the heliospheric shocks is too small to 
map them into x-rays. These groups of waves we will briefly mention below; they will however play a more important 
role in Part 2 of this text when discussing measurements and observations of the various types of shock waves that 
are encountered in the heliosphcrc. 

1. Ion-beam driven instabilities - u) < uo C i 

The shock waves in the heliosphere are magnetised. As long as we are interested in their formation and properties 
we can restrict to low frequency electromagnetic waves in warm plasma. Such waves are excited by plasma streams 
or kinetic anisotropics in one or the other way. The simplest instability known which distorts the magnetic field by 
exciting Alfven waves that are propagating along the magnetic field is the firehose mode. 

a. Firehose mode. The firehose mode is the result of a pressure (or temperature) anisotropy in plasma with the 
parallel pressure Py exceeding the perpendicular P± and magnetic B 2 /2fiQ pressures. |46j gave a simple intuitive 
explanation of this instability based on the insight that the parallel thermal motion of the adiabatic magnetised ions 
along the magnetic field exerts a centrifugal force on the field lines. When this force exceeds the restoring forces 
of the magnetic pressure and perpendicular plasma pressure, the centrifugal force wins and a small excursion of the 
magnetic field starts growing and propagates as a wave along the magnetic flux tube like on a string. The condition 
for instability is 

P|, -P ± > B 2 / M o (74) 

Since the pressure anisotropy on the left means that there is an excess in parallel energy in the plasma, the plasma 
possesses free energy which by the instability is fed into the excitation of Alfven waves with frequencies ll>a *C t-u cil 
transported away with Alfven speed and ultimately dissipated in some way - as expected. The waves excited are 
ordinary Alfven waves, however, and not suited for shock formation. Below we will once more encounter this mode 
in discussing ion beam instabilities. 

b. Kinetic Alfven waves. Excitation of kinetic Alfven waves requires (3 < 1 and a different process. In the 
solar wind the /3-condition is barely satisfied except possibly in the very strong coronal magnetic fields or locally 
(possibly in Corotating Interaction Region boundaries when the magnetic field may become compressed without just 
forming a shock) . Kinetic Alfven waves possess a finite electric field component parallel to the magnetic field which ca 
accelerate electrons. However, the inverse mechanism is also possible that electrons moving along the magnetic field in 
the opposite direction, become retarded by this field component and feed their energy into the kinetic Alfven wave. A 
process similar to this has been suggested by |25j in different context for bouncing electrons in a locally inhomogencous 
magnetic field represented as B(z) = Bq(1 + az 2 ). The electron beam conserves the magnetic moment when moving 
along the magnetic field, interacting adiabatically with the parallel wave electric field for long wavelength 

h 2 r 2 

t? _ _L ci 

" l + *ir£ 

The condition that the electrons form a beam is that at the resonance with the wave the derivative of the electron 
distribution with respect to the resonant parallel electron energy of motion en = m e uj 2 /2b is positive, dF e /de\\ > 0. 
Here the square of the bounce frequency is u> 2 — 2fiuj ce /m e , and /i = m e V 2 /2B is the magnetic moment. The one 
can calculate the growth rate as 

7 kaw (k ± \ e )- 2 (io 4 /bkfv 3 e ), u = k l{ V A (76) 

We see that an electron beam moving in an inhomogeneous magnetic field can excite kinetic Alfven waves. 

Probably more important than this is, however, the interaction of ions which are reflected from a solitary pulse and 
move back upstream ahead of the pulse, as had been suggested by [IB]. The reflected ions will represent a beam that 
is moving against the initial plasma inflow which by itself is another ion beam neutralised by the comoving electrons. 
This configuration leads to a ion beam-ion beam interaction and should cause an instability because free energy is 
present in the two counter-streaming beams. The various instabilities this process may cause have been reviewed by 

c. Kinetic growth rate. Before coming to discussing the relevant instabilities we should briefly mention the waves 
which can be driven by them. We already noted that in thermal plasma most waves will rest in thermal fluctuations. 
Once a wave which is an eigenmode of the plasma is injected it will experience thermal damping until it disappears 
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in the background fluctuations. Hence, a wave that is assumed to grow must overcome this damping which for 
propagation parallel to the average magnetic field Bq in a uniform plasma is given by 
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where the index s identifies the species, v s — y / 2T s /m s is the thermal speed of the species s, V s its average parallel 
bulk drift velocity, lo s = oj ps its plasma frequency, and we have dropped the index r (for real) at oj which in this 
expression is understood as real anyway. Note that we are going to take into account several different ion species and 
thus need an extra index to distinguish between them all. The simplified cold dispersion relation is 
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which determines the approximate real frequency. However, when thermal effect are included, then there is no way 
to avoid the numerical solution of the full kinetic dispersion relation. In the following we will generally refer to such 
numerical solutions. 

It is obvious from the expression for the damping rate, that for sufficiently large average drift velocities V s of species 
s larger than the phase velocity uj/k\\ of the wave this particular species contributes a positive term to the damping 
rate 7 which, when large enough can dominate the entire damping rate. In this way streaming is one way to cause 
instability. In the absence of streaming 7 is independent of V s and is negative for a thermally isotropic plasma. In 
the presence of a temperature anisotropy, however, this may change as we have seen for the firehose mode. The above 
damping rate (77 1 does not account for thermal anisotropy which is, however, simple matter to include [cf., e.g., 121] . 

We note finally that 7 is a resonant damping/growth rate yielding resonant instability driven by small groups of 
resonant particles. The firehose mode is non-resonant since all particles contribute to it. Generally most non-resonant 
instabilities can only be found by solving the full dispersion relation numerically. 

At low frequencies it suffices for our purposes of understanding shock physics to deal with a three-component 
plasma consisting of two ion species and one neutralising electron component which we assume to follow a Maxwcllian 
(thermal) velocity distribution. Moreover, we assume that the drifting ion components are Maxwellians as well. In 
conformity with the above remarks on a resonant instability we assume that the dominant ion component has large 
density 3> Nj,, and the second component represents a weak fast beam of density iVf, propagating on the ion- 
electron background with velocity V> ^ V, ~ 0. Following 21J it is convenient to distinguish the three regimes: cool 
beams (0 < Vb < Vb), warm beams (vb ~ Vb), and hot beams (vb 3> Vb). Figure 12 shows the beam configurations for 



these three cases and the location of the wave resonances respectively the position of the unstable frequencies. 

d. Cool ion beam: Right-hand instability. Assume that the ion beam is thermally isotropic and cool in the above 
sense, i.e. its velocity relative to the bulk plasma is faster than its thermal speed. In this case a right-handed resonant 
instability occurs. In the absence of a beam Vb — the parallel propagating mode is a right-circularly polarised 
magnetosonic wave propagating on the lowest frequency whistler dispersion branch with u) « A:||Va- In presence of 
a drift this wave becomes unstable, and the fastest growing frequency is at frequency lo ~ fcyVJ, — uj C i. This mode 
propagates parallel to the beam, because u> > 0, fcii > 0, and Vb > 0. The numerical solution of this instability for 
densities 0.01 < Nb/Ni < 0.10 at the wave-number k\\ of fastest growth rate identifies a growth rate of the order of 
the wave frequency 7 ~ lu and 



j m ~ u ci (N b /2N e )3 



(79) 



for the maximum growth rate j m , where N e = Ni + Nb is the total density from quasi-neutrality. This instability 
drives waves propagating together with the beam in the direction of the ion beam on the plasma background which 
has been assumed at rest. If applied, for instance, to shock reflected ions then for 2% reflected ions the maximum 
growth rate is ~ 0.2oj ci , and Vb ~ 1.2u> c i/ku , and fen ~ 0.2u — cI/Va which gives Vb ~ QVa- In the solar wind the Alfven 
velocity is about Va ~ 30km/s. Hence the velocity difference between shock reflected ions and solar wind along the 
magnetic field should be roughly ~ 180km/s. The thermal velocity of the ion beam must thus be substantially less 
than this value, corresponding to a thermal beam energy less than Tb <SC 100 eV which in the solar wind, for instance, 
is satisfied near the tangential field line. The solar wind travels at 300-1200 km/s. Complete reflection should produce 
difference speeds twice these values. The above value is thus not unreasonable for travelling shocks, but for bow shock 
reflected ions applies to the quasi-perpendicular portion of the bow shock only. We may thus conclude that this wave 
mode could be excited in the solar wind by shock reflected ion beams near quasi-perpendicular shocks. 

e. Warm ion beam: Left-hand instability. The above instability is present when the ion beam is rather cold. 
When the temperature of the ion beam increases and the background ions remain to be cold, then beam ions appear 
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FIG. 12: The three cases of ion beam - plasma interaction and the location of the unstable frequencies. Shown is the parallel (reduced) 
distribution function i^(fc|itj||), where for simplicity the (constant) parallel wave-number fen has been included into the argument. Right 
handed resonant modes (RH) are excited by a cool not too fast beam. When the beam is too fast the interaction becomes nonresonant. 
When the beam is hot, a resonant left hand mode (LH) i excited. In addition the effect of temperature anisotropy is shown when a 

plateau forms on the distribution function [after [21] . 



on the negative velocity side of the bulk ion distribution and go into resonance their with the left-hand polarised 
ion-Alfven wave. Their maximum growth rate is a fraction of the growth rate of the right-hand low frequency whistler 
mode. Nevertheless it can excite the Alfven-ion cyclotron wave which also propagates parallel to the beam. For this 
instability the beam velocity must exceed the Alfven speed V& > Va- 

At oblique propagation both the right and left hand instabilities have smaller growth rates. But interestingly, it 
has been shown |24| that the fastest growing modes then appear for oblique k and harmonics of the ion cyclotron 
frequency u> ~ nu> c i, with n = 1,2, ... . 

/. Nonresonant ion instability: Firehose mode. When the ion beam is fast and cold it does not go into resonance 
because its velocity is to high. In this case all ions participate in a nonresonant instability which in fact is a thermal 
firehose mode where the ion beam has sufficient energy to shake the field line. This mode propagates antiparallel to 
the ion beam, has small phase speed and negative helicity. This mode has large growth rate for large Nb/N e and 
V^/Va simply because then there are many beam ions and the centrifugal force is large while the beam velocity lies 
outside any resonant wave speed. It is trivial that this instability becomes stronger when the ion beam is composed 
of heavier ions as the larger mass of these increases the centrifugal force effect. 



2. Electron instabilities and radiation - lu ~ Lu pe 

Other than ion beam excited instabilities electron-beam instabilities are not involved in direct shock formation (unless 
the electron beams are highly relativistic which in the entire heliosphere is not the case; we do not investigate relativistic 
shocks in this text anyway). The reason is that the frequencies of electron instabilities are high. However, just because 
of this reason they are crucial in anomalous transport being responsible for anomalous collision frequencies and high 
frequency field fluctuations. The reason is that the high frequency waves lead to energy loss of the electrons retarding 
them while for the heavier ions they represent a fluctuating background scattering them. In this way high frequency 
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FIG. 13: The regions of instability of the electron beam excited whistler mode in density and beam velocity space for two different f} 
compared to the ion acoustic and electron beam modes. Instability is above the curves. The whistler instability has the lowest threshold 

in this parameter range [after [21]. 

waves may contribute to the basic dissipation in shocks even though this dissipation for supercritical shocks will not 
be sufficient to maintain a collisionless shock or even to create a shock under collisionlcss conditions. This is also 
easy to understand intuitively, because the waves need time to be created and to reach a substantial amplitude. 
This time in a fast stream is much shorter than the time the stream needs to cross the shock. So waves will not 
be accumulated there; rather the fast stream will have convected them downstream long before they have reached 
substantial amplitudes for becoming important in scattering. 

When we are going to discuss electromagnetic waves which can be excited by electrons we also must keep in 
mind that such waves can propagate only when there is an electromagnetic dispersion branch in the plasma under 
consideration. These electromagnetic branches in w,k-space are located at frequencies below the electron cyclotron 
frequency oj ce . The corresponding branch is the whistler mode branch. Electrons will (under conditions prevailing 
at shocks) not be able to excite electromagnetic modes at higher frequencies than uj ce . We have seen before that ion 
beams have been able to excite whistlers at low frequencies but above the ion-cyclotron frequency. This was possible 
only because the whistler branch exists due to the presence of electron (because of quasi-neutrality) as a channel for 
wave propagation. In fact, ion cyclotron waves can of this reason also have higher electromagnetic harmonics. For 
electrons under conditions in the shock environment this is generally not possible. Electromagnetic waves excited 
by them propagate on the whistler branch or its low frequency Alfvenic extension [92 . Electron beams thus excite 
electromagnetic whistlers and right-handed Alfven waves. They also excite a variety of electrostatic emissions which 
we will mention later as well. 

a. Whistlers. Whistlers can be driven in two ways, either by an electron temperature anisotropy [29 , or by 
electron beams (or heat fluxes) |21j . In the former case the condition for instability are that the perpendicular 
electron temperature T e j_ > T e u exceeds the parallel electron temperature, and that the parallel energy of the resonant 
electrons £|| = ^rn e V^ > B 2 /2/ioN e exceeds the magnetic energy per electron. 

|21j has investigated the more relevant case of whistler excitation by an electron beam. He finds from numerical 
solution of the full dispersion relation including an electron beam in parallel motion that with increasing beam velocity 
Vb the real frequency of the unstable whistler decreases, i.e. the unstably excited whistler shift to lower frequencies on 
the whistler branch while remaining in the whistler range Ud < u> < lu cc . Both, the background electrons and beam 
electrons contribute resonantly. The most important finding is that the whistler mode for sufficiently large Pi ~ 1 
(which means low magnetic field), Nb/N e and Tf,/T e has a the lowest beam velocity threshold when compared with 
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the electrostatic electron beam instabilities as shown in Figure 13 This finding implies that in a relatively high-/? 



plasma a moderately dense electron beam will first excite whistler waves. In the shock environment the conditions 
for excitation of whistlers should thus be favourable whenever an electron beam propagates across the plasma along 
the relatively weak magnetic field. The electrons in resonance satisfy Vii = (u> — u ce )/kii and, because u> <C u> ce 
the resonant electrons move in the direction opposite to the beam. Enhancing the beam temperature increases the 
number of resonant electrons thus feeding the instability. On the other hand, increasing the beam speed shifts the 
particles out of resonance and decreases the instability. Hence for a given beam temperature the whistler instability 
has a maximum growth rate a few times the ion cyclotron frequency. 



E. "Transport ratios" 

Measured wave spectra are complex and opaque, i.e. it is very difficult from an inspection of their shape to identify the 
wave modes that are present in the plasma volume under investigation. In some rare clean cases one can conclude from 
the observation of a particular maximum in the wave power or the observation of only one single field component which 
wave has been detected. In the general case of broad spectra or mixed spectral maxima and various field components 
lacking clear dominance of one field component it becomes nearly impossible to decide about the waves. Clearly, when 
the spectrum is shapeless power-law one in most cases is dealing with well developed turbulence in which case it 
makes no sense to distinguish and search for single modes. Then one must seek shelter among the well developed 
methods of turbulence. On the other hand, if the spectra indicate the presence of single waves, one would like to 
have some quantities at hand which help identifying which modes one is dealing with. It would be helpful if one 
could measure simultaneously both, the wave frequency and the wave number spectra. This is possible, however, only 
with sophisticated multi-spacecraft constellations. And even then only in the rarest cases the determination of the 
dispersion relation from experiment will be possible. In application of these theoretical arguments to real observations 
one therefore has defined some quantities, called "transport ratios" , which have turned out to be quite valuable in 
helping identity some of the wave modes. Such transport ratios for electromagnetic waves have been given by |21j . 
Below we list the most interesting for our purposes. 

b. Polarisation. The polarisation of a wave magnetic field with respect to wave number k is given by 

P = i°s/b A (80) 

where &s, b A are the components of the magnetic fluctuation field b in the directions S, A of magnetosonic and Alfven 
waves, respectively, i.e. the vector A = k x Bo is perpendicular to the wave vector and the ambient magnetic field, 
while the vector S is perpendicular to k (because of the vanishing divergence V • b — 0) in the plane (k, B ). The 
waves are more magnetosonic or more Alfvenic whether \P\ > 1 or \P\ < 1, rspectively. For ReP > (< 0) the waves 
are right-hand (left-hand) polarised. 

c. Compression. The magnetic compression of the wave measures the relative variation in the parallel magnetic 
fluctuation field 

C B = (6|)/(|b| 2 > (81) 

where the fluctuations are taken at a given pair (w,k). The angular brackets (ab) mean taking the real part of the 
correlation function of the two bracketed quantities. 

d. Parallel compressibility. This ratio together with the compression ratio provides a tool for estimating how 
compressive a wave is. It is defined for species s as 

flg A N s ) 

C||s - (&|) (82) 

e. Non-coplanarity ratio. This ratio measures the fluctuating field component out of the plane (k, B , and is 
given by 

C c =(|^| 2 )/(|b| 2 ) (83) 



/. Alfven ratio. Defining AV^ = h/y/norriiN, where N is the total plasma density, the Alfven ratio is defined as 

R As = (|AV S | 2 )/(|AV A | 2 ) (84) 



Here AV S is the flow velocity of species s. An Alfven wave has AVj = ±V A h/'B , and its own Alfven ratio is R A i = 1. 
The Alfven ratio thus measures the fraction of Alfven waves contained in the near-zero frequency fluctuations. 
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g. Cross-helicity. Hclicity of a wave is another identifier of the wave mode, it is in particular useful for determining 
the direction of propagation of the wave by considering its sign. One has therefore defined a ratio which provides a 
measure of it: 

^ cs = (b.AV s )/[(|AV s | 2 ) + (|b| 2 )] (85) 

Parallel propagating Alfvcn waves have hclicity H ci = — 1. Fast magnetosonic nearly parallel propagating modes have 
also \H ci \ ~ 1, and H ci = for perpendicular propagation. 



V. ANOMALOUS TRANSPORT 



It has been mentioned several time that the evolution of shocks requires the generation of some kind of dissipation. 
Under the conditions of non-collisionality the generation of dissipation must be intrinsic to the plasma. In fact, 
assume that a source of free energy is switched on in the plasma. It is then quite natural to imagine that this available 
free energy will act on the plasma in a way to dissipate itself, distribute itself all over the plasma and ultimately 
transform itself into heat, entropy and create a new thermal equilibrium. Seen from this point of view the occurrence 
of instability is the first step in this chain of processes directed toward thermal equilibrium. Putting an obstacle into a 
fast but otherwise thcrmalizcd plasma stream is clearly a way of providing free energy, because seen from the frame of 
the obstacle the plasma is not in equilibrium; there is a high velocity difference between plasma and obstacle and thus 
a large amount of free latent energy available in the system. Dissipation, however, requires small scale interactions 
between the plasma constituents. Macroscopic motions with their long scales contribute to large scale structures but 
do not directly act on the microscopic scales. In order to affect the particle motion and contribute to friction among 
the particles small scale processes have to be called for. 

These processes take mostly place on the electron scales. Moreover, since it is much easier to excite fluctuations 
in the electric field than in the magnetic field, these interactions arc electrostatic. In the following we consider the 
electrostatic fluctuations which are expected to contribute to the generation of anomalous dissipation. These processes 
can be divided into those which are not affected by the presence of an external magnetic field called unmagnetized, and 
those where the external magnetic field must be taken into account in the particle motion, i.e. magnetized processes. 
The distinction is made by the relation between plasma and cyclotron frequencies. 

When Va <C C\Jm e jm^ the electron cyclotron frequency is much less than the electron plasma frequency, uj pe » uj ce 
and one is dealing with an unmagnetized case. Also, when the entire dynamics is restricted to the direction parallel 
to the magnetic field, the problem can be considered to be unmagnetized. The complete unmagnetized dispersion 
relation including all species s and their drifts V s is 

i + X>.Kt) = o, A-,(., k )^A|^M m 

The function K s is the susceptibility contribution of species s with average distribution function F 0s . In a Maxwellian 
component plasma K s = — (l/2k 2 X 2 - )s )Z'(^ s ) can be expressed through the plasma dispersion function Z(£ s ), with 
\d.s the Debye-length, ( s — (uj — k • V s )/^/2kv s , and thermal speed v s of component s. Note that here the sign of 
charge is included in the cyclotron frequency, i.e. for electrons — uu ce , for ions +oj ci . 

When the plasma is magnetized, which applies to all other cases, the susceptibility becomes more involved. For 
Maxwellian components it reads 



v2|K|||t> a 



Here rj s = (k±r cs ) 2 with Ii(r] s ) the order-/ Bessel function of imaginary argument. 

This last expression suggests that magnetized electrons will support electron cyclotron harmonics lu) ce , which, 
when purely perpendicular, are Bernstein modes. There are also ion Bernstein modes luj ci , but (with the exception 
of injection of a localized perpendicular ion beam) they do not play any susceptible role in transport. Exciting both 
electron and ion Bernstein modes requires beams perpendicular to the magnetic field, in which case when applied to 
shocks the lower hybrid instability will become the most important agent in generating anomalous dissipation. 

As with the electromagnetic waves discussed in the previous sections the electrostatic dispersion relations can be 
solved analytically in closed form only in particular simple cases, and one has to retreat to a numerical approach. In 
the following we present the few most important cases and their effect on the generation of anomalous dissipation 
keeping in mind, however, that most collisionless shocks in the heliosphere are supercritical, and anomalous dissipation 
does not contribute substantially to their evolution and maintenance. Dealing in Chapter 5 with supercritical shocks, 
particle reflection is the dominant dissipation process, and we will refer to the instabilities of the previous section. 
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A. Electrostatic wave particle interactions 



The wave friction term Eq. ( 26 1 requires the determination of the fluctuation amplitudes of the wave modes that are 
responsible for causing anomalous friction. Again, the first step is to identify the unstable wave modes. In the second 
step we will then either have to determine the saturation level amplitudes of these waves or to consider their further 
interaction with particles or other waves. 

Here we list instabilities of interest in the anomalous dissipation process only. When discussing application to 
observations we will later in passing also mention instabilities which are involved in radiation from shocks. However, 
radiation provides no substantial energy loss, and the dissipation caused by radiation under the conditions of the 
nonrelativistic shocks in the hcliosphcrc is completely negligible and docs neither affect shock formation nor shock 
structure. One possible exception are shocks in the solar atmosphere which sometimes are accompanied by x-ray 
emission which, however, is not a genuine unstable plasma process in this case. Mentioning of such processes in 
relation to solar coronal observations is deferred to Chapter 8. 



1. Unmagnetized electron and ion instabilities 

The unmagnetized dispersion relation in a Maxwellian component plasma consisting of one electron and one ion 
component has two solution, electron plasma waves or Langmuir waves to 2 — to 2 e + 3k 2 v 2 at k\\\De <C 1 and ion- 
acoustic waves. The Langmuir wave can be driven unstable by a parallel electron beam of velocity V e > 3v e . 

a. Ion-acoustic waves. More interesting are ion-acoustic waves. In the absence of any difference velocity between 
electrons and ions these are strongly damped plasma waves propagating along the magnetic field with dispersion 
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These waves become dispersive at larger fcy and for large electron temperatures with concave dispersion curve. At 
small parallel wave numbers their dispersion is linear resembling sound waves with velocity c 2 ia ~ [T e + 37i)/raj. 
Obviously for dispersion to compensate nonlinearity large electron temperatures T e > 3T^ are required in which 
case also the damping is small. Ion acoustic waves if managing to overcome damping are therefore a candidate for 
electrostatic unmagnetized shock formation: Their dispersion favours shock ramps with slow moving wave trails and, 
as we will see later, they also can contribute to dissipation thus satisfying all conditions as candidates for shock 
formation. However, such shocks are purely electrostatic and do not affect the magnetic field. They will therefore 
only be of interest in sub-structuring magnetized shocks, possibly contributing to the formation of subshocks of short 
wavelength of the order of several Debye lengths. 

In fact going from the kinetic description to the fluid description and writing down the continuity and momentum 
conservation equations for ion acoustic waves for parallel propagation and one-dimensionality and combining it with 
Poisson's equation 
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under the well justified assumption that for the low ion-acoustic frequencies the electrons behave as thermalized hot 
Boltzmannians with density N e — N exp(e0/T e ) depending exponentially on the electrostatic potential </>, assuming 
quasi-ncutrality and localized stationary solutions, we manipulate all these equations into the Sagdeev pseudo-potential 
form 
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As usual the subscript 1 refers to values far upstream of the localized solution. In the absence of dissipation this solution 
for S < yields solitons of maximum potential amplitude <f> m found from setting the bracket to zero. (f> m corresponds 
to a maximum compressive amplitude N rn = N exp(ecj) m /T e ). With ion acoustic Mach number M. ia = V\/ci a it is 
found that solitons exist only in the supersonic regime M.i a > 1. The soliton speed can be expressed through the 
maximum potential (or density via Boltzmann's expression) as 
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FIG. 14: Left: The threshold drift speed for the electron current driven ion-acoustic instability as function of the electron to ion 
temperature ratio. For comparison the threshold for the parallel electron current driven ion cyclotron instability in a magnetized plasma 
is show for a dense plasma with LU pe /uj ce = 10 [after l30l . Right: Buneman electron drift-ion two-stream instability, frequency and growth 
rate as function of wave number at large electron current drifts when the ion-acoustic instability has changed to the Buneman mode. 



The condition on the maximum compression amplitude N m in this expression simply requires that the potential 
energy must be less than the initial flow energy. This sets a limit on the possible Mach numbers A4% a > 2e</> m /T e or, 
when combined with the definition of the latter, A4f a /2 > ln(N m /N). 

Since the denominator in the former equation must be real, this condition requires that 

HN m /N) + 1 - N m /N < 

Expressing herein the density ratio through the Mach number, one finds that for ion-acoustic solitons to exist the 
Mach number is limited to values below a surprisingly small critical Mach number 

M ia < Mf* ~ 1.6 

Ion acoustic solitons cannot exists for Mach numbers exceeding -M?"*. For higher inflow velocities V\ the ion acoustic 
soliton will cither not evolve or break down. 

The range of possible Mach numbers is rather limited which is simply due to the fact that for higher speeds 
the dispersion is unable to sustain a stationary state. Any solution will be non-stationary, wave like or unstable. 
Concerning the formation of shocks one in addition to dispersion requires that dissipation is produced. Since it is 
known that ion acoustic waves are Landau damped with damping rate 7l,ioj one can argue that Landau damping will 
cause a shock profile on the ion acoustic soliton with downstream state different from the upstream state [I3J |3S] • The 
downstream density is then found to be iVjj ~ N\ exp[2(Mi a — 1)] and will exhibit trailing oscillations, as has been 
discussed above. However, Landau damping takes time, and therefore the general argument applies to this kind of 
shock formation that the damping will not have time to work for large Mach numbers. Thus the damping argument 
applies only to subcritical shocks of Mach numbers smaller than M.i a ~ 1.6. Such weak (electrostatic) ion acoustic 
shocks can indeed evolve and may contribute to sub-structuring of stronger supercritical shocks in the region where 
the Mach number has already dropped to values below the critical. 

[4"6] has favoured reflection of inflowing particles from the leading edge of the soliton [35] over Landau damping. This 
reflection affects ions with energy less than the soliton potential <fi m and causes oscillations of long wave length. More 
important is that the reflected ions form an ion-ion beam configuration and are thus subject to the ion instabilities 
discussed previously yielding waves which may generate dissipation but do also propagate upstream of the shock 
where they cause wave particle interactions and retard the inflow ahead of the shock. 

b. Electron current driven ion acoustic instability. So far we have not asked for the reason of an ion acoustic 
wave to grow. This can be achieved in the simplest way by letting one of the plasma components drift with respect 
to the other. If in a two-component plasma the electron drift with respect to the ions they effectively carry a current 
j|l = — eiWeH which is in most cases - but not necessarily - along the magnetic field. Here we assumed again quasi- 
neutrality N e — Ni = N which dispenses us from considering space charges and solving Poisson's equation. In this 
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case assuming weak growth such that we can apply the general instability theory with 7 <C u) the growth rate of the 
ion acoustic wave Eq. ( 88 1 becomes 
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Instability sets on for V e > u)/k\\ when the electron velocity, which is the current drift velocity, exceeds the phase 
velocity of the ion acoustic wave, i.e. when - approximately - V e > Cia. These waves have relatively long wavelength 
&||A_D e -C 1. The threshold for marginal stability of these electron current driven ion acoustic waves can be obtained 
from setting "/i a = 0, yielding 
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The second term in the brackets results from the Landau damping of the ion acoustic waves. It is seen that this 
term disappears for hot electrons with T e 3> Ti thus lowering the threshold for instability to its marginally smallest 
value V e — Cia- The threshold is shown graphically as function of the temperature ratio in Figure [T4) The threshold 
is measured in ion thermal speed Vi and is quite high. Moreover, the electron temperature must be high implying 
that the electron distribution must be hotter than the ion distribution. In the solar wind this is usually satisfied. 
The physical reason for the electron temperature to be high is that the distribution must have a positive slope in v in 



the region of overlap with the cold ion distribution for resonant instability, dFo e /dv 
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> 0. There must be more 



fast than slow electrons in the phase velocity frame of the wave in order to push the wave to higher momentum and 
energy, i.e. causing instability. 

It is clear that this kind of interaction between the ion acoustic wave goes on the expense of the motional energy 
of the resonant electrons. Hence one expects that ion acoustic waves retard and scatter the current electrons thereby 
reducing the current flow. This resembles collisional friction which the resonant electrons experience and can thus 
be interpreted as the generation of an anomalous resistance in the plasma. An interpretation like this has been put 
forward by |46j and has been elaborated in depth afterwards. Below we will return to this theory. 

c. Buneman electron current two stream instability. Returning to Figure |T4"| we observe a change in the threshold 
curve for V e /vi > \J mi/m e T e /Ti. Here the electrons become thermally slow with respect to their drift velocity, 
V e < V e . The ion-acoustic instability under this condition changes into the Buneman instability which is an electron 
current fluid instability and is also known under the name electron current drift or electron-ion two stream instability. 
It has been discovered by [5] and favoured for application in shock physics already by It should be noted that 
the transition from ion-acoustic to two-stream instability has been investigated in depth by |15j . 

Now, treating the electrons as cold the kinetic effects disappear, and the complex dispersion relation of the Buneman 
instability becomes 
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Note that under these conditions the weak growth theory cannot be applied anymore. Instead one must find the 
growth rate from the complex solutions of this quartic expression. Fortunately, this equation can be solved since for 
resonant electrons the third term becomes dominant. The instability has real frequency lu ~ fcj|V e and maximum 
growth rate of the order of the ion plasma frequency 
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Figure[l4] on its right shows the Buneman frequency and growth rate for a case of very large electron current drift 
V e = 600 Vi in dependence on the wave number fcy Xdz- The maximum growth of the instability is close to kuX^i ~ 0-9 
at short wavelengths slightly larger than the Debye length. The growth rate of this instability is very large. This 
implies that the instability is very strong and grows very fast thereby consuming a substantial fraction of the current 
streaming energy. The nonlinear treatment of this instability is of particular interest for shock physics. Since the 
instability grows so fast it makes little sense to treat it analytically for reasons which will become clear when dealing 
with the application of numerical simulations to shocks. The physics involved into the two stream instability can 
be described as follows. Both the electrons carrying the current and the ions are cold. The situation is thus two 
stream, and the instability is not resonant but reactive with all particles participating. This is the reason for its 
strength and rapidity. Because it consumes a fraction of the bulk flow energy of the electron current, the current 
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FIG. 15: Left: Regions of existence of electron and ion acoustic modes in the density/temperature plane. Here N = N c + Nf, = Ni is 

the total quasi-neutral density. The temperatures given as the hot to cold temperature ratio. The electron acoustic mode extends to 
larger cold density and higher hot temperatures than the ion acoustic mode. Right: Unstable versus stable domains for electron and ion 
acoustic modes in the density /drift velocity plane for the case when the electron acoustic mode is destabilised by a drifting hot plasma 

carrying a current. 



becomes decelerated, and the energy is going mainly into the electrons which are heated by the instability until the 
instability stabilises when V e < v e . Then the ion-acoustic instability takes over. The Buncman two stream instability 
is thus accompanied by a burst in electron temperature and a rapid decrease in current. However, the final state of 
the instability is not a stationary state because the ion acoustic instability continues to grow with its own dynamics, 
possibly ending up in the formation of solitons when the Mach number has sufficiently decreased by current relaxation 
and heating, the former decreasing V e , the latter increasing Cj a . But even during the blow-up phase of the two-stream 
instability structuring similar to soliton formation occurs. This can only be inferred from numerical simulation. 

d. Modified two-stream instability. This is a variant of the two-stream instability driven by a relative drift between 
electrons and ion across the ambient magnetic field B [37 . The dispersion relation for the modified two-stream 
instability is 
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This expression is written here in the electron frame of reference and with angle Qbu- One recognizes that the first 
term is the ordinary Buncman two-stream term. However, for oblique propagation 8g n ^ 90° and 0_b„ =/= 0° the 
two-stream mode couples to the whistler mode. It is this coupling which makes the modified two-stream instability 
interesting for shocks. Dispersion curves and growth rates are shown in Figure 16 for Qsn — 60°, and Vi = Va, and 
an artificial mass ratio m,ijm e = 80 which has been taken in view of numerical simulations to be discussed later. 

The modified two-stream instability (MTSI) operates also for relative drifts smaller than the electron thermal 
but larger than the ion-acoustic velocity and even for Ti ~ T e which makes it potentially important if only such 
perpendicular drifts can be generated. The unstable frequency is in the range of the lower-hybrid frequency. Hence 
the ions can be taken unmagnetized with strongly magnetized electrons. However, it requires oblique relative electron 
drifts since for perpendicular drift the instability disappears meaning that the unmagnetized ions propagate under an 
angle to the magnetic field while the magnetized electrons move only parallel to the magnetic field. 

e. Electron-beam electron- acoustic instability. A mode very similar to the ion-acoustic mode is the electron acous- 
tic mode. Here the role of the ions is played by the cold (subscript c) electron background. In presence of another hot 
(subscript h) electron component the physics becomes very similar to the ion-acoustic wave, and a weakly damped 
resonant mode with real dispersion relation 
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can propagate in the range of wave numbers indicated on the right. This weak damping can, like in the case of ion 
acoustic instability, overcome if the hot electron component drifts or if another electron beam is injected into the 
plasma. In the former case with ions and cold electrons at rest the drift implies current flow. The regions of existence 
and stability of the electron acoustic instability are shown in Figure[l5]in comparison with the ion acoustic instability. 
Growth rates have been found numerically to be quite large, the order of the ion plasma frequency and thus similar 
to the Buneman mode. On the other hand, numerical experiments have not detected strong electron acoustic waves. 

Inspection of the electron acoustic dispersion relation shows that the mode is dispersive as well. For very low cold 
electron temperatures this dispersion is similar to ion acoustic waves suggesting that electron acoustic solitons could 
be formed in a similar way. In fact, such solitons have been calculated analytically for different parameter regimes 
[14] . However, observations do not seem to support their existence. Also numerical simulations |36j have not shown 
the formation of such solitons. Clearly, electron acoustic waves can have a different dynamics because of the fast 
reaction times of the electrons, leading to rapid thermalization of the plasma. These questions still remain to be open. 
In any case, if electron acoustic waves exist, the electron acoustic mode is quite well suited for plasma heating in 
shocks where plasmas of different temperatures mix. Its advantage is also that it proceeds on a very fast time scale 
close to the inverse of the plasma frequency. It is hence well suited for fast production of anomalous dissipation of 
energy. Moreover, since it very obviously damps rapidly it may act as an agent of about immediate transformation 
of excess energy in the electrons into heating electrons. None of these question has been understood nor answered 
properly at the time of writing. 

/. Ion-beam ion-acoustic instability. An instability very similar to electron beam driven electron acoustic insta- 
bility is its ion pendant when a cold ion core and hot ion beam in presence of a neutralising electron background 
become unstable |23j . The mode excited in this case is again the ion acoustic mode, this time propagating at very low 
phase velocities to/k ~ Ci a (Ni_ c — N — i,b)/N e less than Ci a . In this range the waves have no dispersion. This wave 
is, however, very easy to destabilise because of its low beam velocity threshold which can lie even below the thermal 
speed of the ion core component. On the other hand the growth rates of this instability are very low. Measured 
as before for the electron-acoustic, ion-acoustic and Buneman instabilities in terms of the ion plasma frequency, its 
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maximum growth rate is just of the order of a few per cent of 7 m ~ 0.3w p i, very small compared to the growth rate 
of the electron-current driven ion acoustic instability which was of the order of j m ~ uj pl . It should thus be less 
important in anomalous processes. 

g. Electrostatic ion- cyclotron harmonic instabilities. So far we have considered only unmagnetized instabilities. 
We now turn to listing the most important - in view of the anomalous processes we have in mind when considering 
shock formation - magnetized electrostatic instabilities. These instabilities occur when a magnetic field is present in 
the plasma and the electrons must be considered to be magnetized. 

The frequencies of the waves in question will therefore be well below the electron cyclotron frequency u> <C oj C c 
falling into the range of and above the ion cyclotron frequency oj > oj C i. According to linear theory the magnetized 
modes in this range occur at harmonics u>i ~ luj C i of the ion cyclotron frequency. We are thus dealing here with with 
electrostatic ion-cyclotron harmonics. 

Following a suggestion by |13j ion cyclotron instabilities have been proposed as generating anomalous collisions 
by [30 who advocated their importance because of their apparently lower instability threshold than the ion acoustic 
instability for electron current shown in Figure [T4| In fact ion-cyclotron harmonic waves oblique to the magnetic field 
can become unstable in the presence of an electron current flowing along the magnetic field. The reason is that - 
in contrast to the earlier mentioned strictly perpendicular Bernstein mode resonances - such oblique waves posses a 
field-aligned projection component of their electrostatic field which can resonant with the parallel current electrons 
via the Doppler-shifted resonance condition A)||V|| e = ui ± luj C i of which only Landau resonance I = is effective for 
\oj — fci|V^|i| < v2|&|||«e- The ions contribute only a weak resonant damping of the waves for I = — 1. Note that the 
obliqueness of propagation of these modes in contrast to bernstein modes implies that the resonance does not occur 
precisely at the harmonics but lies somewhere in between in the harmonic bands where the dispersion curves have 
particular geometrical forms [cf., any book on basic plasma wave physics, e.g.,[2]. Its precise location depends heavily 
on the exact prescribed conditions of the plasma and current velocity V e t\, and no analytic expression can be provided. 

Strongest weak instability occurs in the harmonic range 1 < I < 2 with growth rate 7 oj and k\\r C i ~ 1, i.e. 
wavelengths comparable to the thermal ion gyro-radius in the frame of the ions and propagation angles 6 < 85°. 
The velocity threshold decreases with increasing electron temperature T e /Ti simply because more electrons go into 
resonance then. But for large ratios it is taken over by the ion acoustic instability as seen from Figure[l4] (Note that 
in the solar wind/bow shock system, for instance, the ratio is about T e /Ti w 10 changing across the bow shock to 
become T e Ti < 1; hence one may expect the first electrostatic ion-cyclotron harmonic to be present wherever parallel 
electron currents flow in the shock on the upstream side, while downstream neither current driven ion cyclotron nor 
ion acoustic instabilities should exist.) 

We have already noted that ion beams can also excite ion cyclotron harmonic waves. Parallel beams excite similar 
waves with similar properties like parallel electron currents for propagation angles 0° <C 9 < 90° rather close to 
perpendicular. Perpendicular ion beams, on the other hand excite a broad spectrum of ion cyclotron harmonics on the 
background ion component depending on beam energy and background temperature. These excitations are restricted 
to a distance of the order of the ion gyro-radius r < r c i,b of the beam ions only, because at larger perpendicular scales 
the beam ions are themselves magnetized. Thus this kind of excitation is of importance merely when heavy ion beams 
penetrate the plasma, for instance when a heavy ion beam component of large perpendicular energy is reflected from 
a supercritical shock wave back upstream, or when ions become highly accelerated in interaction with the shock and 
penetrate across the shock onto the downstream side where they have much larger gyro-radii than background ions. 

h. Electrostatic electron- cyclotron instability. This instability is the pendant to the former ion cyclotron instabil- 
ity at the much higher frequency perpendicular (or oblique) electron cyclotron harmonics (for purely perpendicular 
propagation these are the Bernstein modes) . It is driven by the coupling between a sufficiently fast ion beam and the 
electron cyclotron harmonics at perpendicular wave numbers k± > and can also be driven unstable when reflected 
ions are present for instance in the foot of the supercritical perpendicular shock. 

2. Electrostatic drift instabilities in inhomogeneous plasma. 

The last group of instabilities we will refer to in the context of shock physics are instabilities that are caused in 
presence of plasma inhomogeneity. Three basic kinds of plasma inhomogeneities can be identified: density N(pz), 
magnetic B(x), and temperature T(x) real space inhomogeneities. The dependencies of these (average) quantities on 
space may in most cases not be independent. In the following, for the purposes of investigation of their effects on 
wave excitation, we will take them as being in fact mutually independent. For simplicity we will take into account 
only density inhomogeneities perpendicular to Bo on scales much larger than the wave length. We also assume quasi- 
linearity which is justified because under the assumption of weak gradients the effects of the inhomogeneity will be 
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weak as well and thus cause only slow wave growth, 
perpendicular direction x up to first order and write 

N s (x) = N s0 (l + e N x), 



In this case one can expand the density with respect to the 



where 
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The effect of the inhomogeneity is that the magnetized particle component behaves adiabatically and starts performing 
a diamagnetic drift motion Vn s y = (eat^s/ ± oJcs) y perpendicular to the magnetic field and density gradient into ±y 
direction, depending on the sign of the particle charge. The ±-sign in the denominator indicates that the cyclotron 
frequency is taken here including the sign of the charge. Drift motions of this kind cause a perpendicular drift current 
jd v to flow in the plasma because particles of different charge sign move in opposite directions. The waves excited 
under such conditions satisfy the modified dispersion relation 
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The new susceptibility K s (en) is still local, i.e. it changes on scales larger than the gradient scale, but locally depends 
only on ejy. It is of the same for as Eq. (87) with two differences: because of the occurrence of the finite perpendicular 



drift Vd y y the frequency in the factor in front of the sum in (87 1 is exchanged with the Doppler-shifted frequency 
lo — > lj — k±VMs- In addition a new drift term 
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must be added inside the brackets. Instabilities resulting from this dispersion relation in the weak instability limit 
bear the general name of drift or universal instabilities. They resonant with the drift motion and have frequencies 
u> ~ k±VNi and long wave lengths satisfying k±r C i < 1. for smaller wavelengths these waves are highly dispersive and 
thus can form nonlinear structures. 

This mode becomes particularly interesting and important in shock physics when the drift speed is so large that 
the ions can be considered as unmagnetized. This happens because the frequency of the drift mode increases with 
Vn and quickly exceeds the ion cyclotron frequency. In this case the frequency of maximum growth is close to the 
lower- hybrid frequency u;/, and the drift mode becomes a lower- hybrid drift wave. 

a. Lower-hybrid drift instability. For unmagnetized ions, when the frequency of the drift wave is ui 3> t-u C i, the 
ion cyclotron frequency can be neglected and the susceptibility simplifies 
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The positive slope on the distribution function which is responsible for instability is in this case on the perpendicular 
part which depends on the drift velocity. There the maximum of the distribution is shifted out of the origin to the 
location of the drift velocity. The real frequency and growth rates are given by [22] 



k±V Ne 

7 



1 



Ti ) I(7 h ) 

irT e V Ne e"° 
2 Ti Vi Io(rie) 



- 1 



Te 



(99) 
(100) 



The unstable wave propagates antiparallel to the direction of the electron gradient drift, i.e. in the direction of the 
electric drift current. Maximum growth of this lower hybrid wave has been found at long wavelengths fc_L,m?*ci ~ 
yj mi/m e over a relatively broad frequency range close to the lower hybrid frequency lu a- t — lu ce \J m e /rrii. 

The reason for the lower-hybrid drift instability to just excite the lower-hybrid frequency is that it is a perpendicular 
two-stream instability similar to the Buneman mode that is, however, driven by the bulk velocity difference between 
electron and ion gradient drifts. It thus propagates on the background ion component. In fact in the (u>, fcjj-plane the 
drift-beam mode u>d — k± VVe couples to many ion cyclotron harmonics thereby exciting almost all of them. However, 
largest growth occurs in the harmonic dispersion band that contains the lower hybrid frequency which is a strong 
plasma resonance. 

The lower-hybrid drift instability is the strongest in the family of the electrostatic ion cyclotron instabilities. As a 
two stream instability its maximum growth rate remains still modestly large being of the order of a few ion cyclotron 
frequencies, "fih.m — (1 — 3)w C i. In a plasma leaving sufficient time r n ;, sat 3> w" 1 for growth and saturation it may 
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well play a substantial role in generating dissipation. We will see later that this instability indeed provides the highest 
so far inferred from instabilities anomalous collision frequency which turns out to be of the order of the lower-hybrid 
frequency itself. Its relevance in application to collisionlcss shocks is however questionable, because of the above 
argument. 

In addition, the lower-hybrid drift instability appears to stabilise under fj > 1 conditions |69j . In application to 
shocks this restriction, if it translates into the nonlinear regime, presents a severe barrier to the use of the lower hybrid 
drift instability as generator of anomalous resistance, dissipation and entropy generation. 



B. Anomalous resistivity 

Resistivity is defined via the Drude formula r\ = v /eoUJp e , with v = a c Nv e the collision frequency. The latter, under 
collisionlcss conditions, becomes the anomalous collision frequency v a and is the quantity containing the interaction 
between electrons and the nonlinear wave fluctuations. 

This becomes obvious when realising that the Spitzer-Coulomb collision frequency vq ~ uj pe /N\^ De is proportional 
to the ratio of the plasma wave fluctuation level in thermal equilibrium Wth = ^ e o( e th) t° thermal energy, vq ~ 
UpeWth/NT e . Under saturated instability conditions it is then reasonable to assume that the actual fluctuations (e 2 ) 
replace the thermal fluctuations in this expression which yields the Sagdeev formula 
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The problem is thus reduced to the determination of the nonlinear saturation level of the unstable wave spectrum. Its 
determination requires knowledge of the electric current j ~ —e(NV e ) as a functional of the electric wave fluctuation 
field e. (We are speaking here only of electrons since in collisionless plasma electrons - because of their much faster 
mobility than ions - are the particles which carry the electric current. The electrons feel the friction of the waves 
and become retarded by anomalous collisions thereby dissipating the kinetic energy of the current and contributing 
to collisionlcss Joule heating of the plasma.) The evolution of the electron current is - in principle - given by the 



electronic part of Eq. (24 1, respectively Ohm's law ( 30 1 , if on the right-hand sides the average anomalous electronic 
friction terms (26 1 are added, since these are the crucial terms containing the wave-particle interactions. 

We are interested only in the parallel collision frequency here. In the nonlinear stationary state the time derivatives 
can be neglected. In order to obtain a first expression for the parallel anomalous collision frequency v a we assume 



that the last term in Eq. ( 24 ) is of the form 

- — NVn = I dv 3 vnC e 



m. 



When inserting from Eq. ( 26 ) , keeping only the parallel component and the electric part which in the microscopic 
interactions dominates over the electromagnetic fluctuations since these affect only frequencies below w c j. One obtains 
for the anomalous collision frequency 

*a*-^V||<W aot > (102) 

which is general but still preliminary. It requires knowledge of the average power density of the electric wave field 
{Wsat) — 5 e o(l e2 |) which here can be of arbitrarily large amplitude and arbitrary spatial structure. Remember that 



the only condition implied in deriving ( 26 ) was that the fluctuations were fast both in space and time compared to 
the slow changes in the plasma background quantities. In this sense the parallel gradient operator Vy acts on the 
slow variability of the wave power. 

In order to proceed, another equation is required which determines the evolution of the wave power. This lacking 
equation can only be formulated in Fourier space (w,k) and should contain all the nonlinear interactions and thus 
cannot be of general nature. A simplifying assumption is that it describes the nonlinear evolution according to kind 
of wave-kinetic equation 

-W k ~-^+[Y-^\ ■VW k = 2 7 (Lu,-k,W k )W k + ... (103) 

We assume the system has reached stationarity such that the wave spectrum W k does not evolve with time anymore. 
In this case the left-hand side simplifies, and we can express the spatial derivative of the stationary wave spectrum as 
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which after transforming from Fourier into real space and averaging over the fast fluctuations can be used in the above 
expression (1021 to express the anomalous collision rate through the average wave power. Such expressions will be 
given below. Usually the current-drift speed |Vm| 3> \du)/dki\\ is much larger than the wave group velocity, and the 
latter can be neglected. This yields the inverse square dependence of v a oc | Vj| |~ 2 on the current drift velocity. 

The derivation of an anomalous collision rate contains a number of crucial assumptions of which the most severe 
concern the simplifications in the terms in the kinetic wave equation (1031 where we suppressed a refraction term 
(Vw) • (dWk/dk) - which becomes important in strongly inhomogeneous plasmas like in a shock ramp - and neglected 
all terms on the right with the exception of the growth rate term. Since the dependence of 7 on the wave amplitude has 
been retained, some generality in the nonlinearity is still retained. Stronger simplifications are made when restricting 
to pure quasilinear theory. In this theory the nonlinear dependence of the growth rate on the wave amplitude is 
dropped. This confirms to the conventional approach to anomalous dissipation. The most elaborated quasilinear 
theory (valid for any direction of propagation including electromagnetic contributions) can be found in [60J where 
also the application to one particular mode (the lower-hybrid drift mode discussed below) is given and a rudimentary 
contribution of the neglected Coulomb collision term is retained. 

In any case the mechanism of saturation of the nonlinear wave field must be known in order to obtain a useful 
practical expression for v a . In the following we review only the three wave modes that contribute strongest to 
anomalous resistance, ion-acoustic, Buncman two-stream, and lower- hybrid drift modes. 



1. Nonlinear evolution of waves 

Quasilinear theory (47J [60] is so far the simplest and most effective approach to the calculation of anomalous collision 
frequencies. This approach uses the linear growth rate 7(k) of the instability yielding the simplified formula 
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It requires in addition knowledge of 7(k) and the drift velocity V e = Viiz. 

a. Anomalous ion acoustic collision frequency. For instance, from this expression, assuming Vn > Ci a and 
kmax^De ~ lj and 7i a ~ ujVi\/v e which holds for the ion acoustic instability in the large drift limit, one just obtains 



the above Sagdeev expression ( 101 1 for the anomalous collision rate which is good for application when the saturation 
level of the instability is measured. 

Ion acoustic waves saturate by scattering off thermal plasma ions according to the resonant scattering process 
lj — lo' = (k — k') • Vi where the prime ' indicates the frequency and wave numbers of the scattered wave. This process 
has been used by (46] to explicitly calculate the ion acoustic anomalous collision frequency 
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which holds for large electron current drifts Vji 3> ci a and for narrow angles 9 > 0. Actually, experiments have shown 
that this expression overestimates the anomalous resistance suggesting that anomalous collisions are less effective than 
theory predicts. More precise theories than the above estimate have been developed by |53] , and have been reviewed 
by 07]. 

b. Anomalous two-stream collision frequency. As noted earlier, the two-stream instability is a very strong insta- 
bility causing large current momentum losses and rapid plasma heating. It switches itself of during evolution and 
will therefore not be a stationary instability. It causes, in fact, different effects which probably destroy its direct 
importance in collision generation. 

It has phase speed oj/k ~ v e yjm e /mi substantially below electron thermal speed implying that it stays relatively 
long in the volume of excitation which supports its effect on the local plasma. Its wave energy density, from simple 
consideration is less than ion energy Wt s <Ti. It yields a large theoretical collision frequency 

VaM — Wpi > V a ,ic With V ailc ~ 0.3/c||U e < LU ci (107) 

of the order of the ion plasma frequency and several orders of magnitude larger than Spitzer-Coulomb collision 
frequency. As such it is a serious candidate for generating anomalous dissipation and heating in shocks whenever 
a two-stream situation is encountered. This is indeed frequently the case as we will sec in the supercritical shock 
Chapters 4 and 5. 



If in the above Eq. ( 107 1 we compare the two-stream collision rate, for instance, with the ion-cyclotron collision rate 
v a ,ic that had been favoured by |30| . We find that it has far larger growth rate than the ion cyclotron wave and will 
thus always dominate when the current is strong, Vfi > v e . At weaker currents the ion-acoustic will be faster (because 
of the slowness of the ion-cyclotron instability) as long as the magnetic field is weak and Vji S> Cj . 



41 



c. Anomalous lower-hybrid drift collision rate. As noted earlier, the lower-hybrid drift instability is particularly 
important as it is the exceptional representant of a fast growing (universal) drift wave instability which in the presence 
of gradients in plasma can always be expected to grow. Clearly shocks are a particularly good candidate for such 
an instability because of the steep density and magnetic field gradients occurring in compressive shocks. Moreover, 
in application to shocks, other than at pressure equilibrium boundaries like the magnetopause [see, e.g.. 1521 155] . the 
magnetic gradient adds positively to the growth rate. 

In calculating the quasilinear saturation level of this instability makes use of the wave number at maximum growth 
kmax^rji = 2/(1 + ujp e /uj^ e ) and d(uie)/doj = 1 + ujp e /ui^ e ^ 1 in dense plasma like the vicinity of the bow shock, for 
instance. The saturation wave level [7 then follows after solving the quasilinear diffusion equation to be fraction of 
the drift energy 

W sat ,lh ^ ^T— 2 / 2 ^ ( Wii ) 

From here the anomalous collision rate follows as 
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proportional to the ratio of ion-gyroradius to density gradient scale length Ln = |VlnA^| _1 . This growth rate is large 
for steep gradient scales close to the ion gyro-radius, a condition that holds in the shock ramp. The propagation of 
electrons in the lower-hybrid drift case is perpendicular to the magnetic field since the electrons perform a diamagnctic 
drift which constitutes the electric current. In the above one-dimensional theory everything has been reduced to the 
coordinate parallel to the current. In a shock wave this current will flow along the shock surface, perpendicular to 
the magnetic field and shock ramp density/field gradient. 

The anomalous lower-hybrid drift collision frequency is large and renders the lower-hybrid drift instability a viable 
candidate for generation of anomalous dissipation in shocks if only the condition that maximum growth is found for 
(3 < 1 can be circumvented. Shocks in the heliosphere in most cases satisfy the inverse condition (5 > 1. Whether this 
indeed provides a serious restriction has not yet been clearly verified, however, neither theoretically nor in numerical 
simulation or observation. 

Recently |60j reviewed the theory of anomalous resistivity for the lower-hybrid drift and two-stream instabilities 
and derived some (slightly) more precise (but considerably more involved) expressions than the formulae given above. 
They also included electromagnetic effects and arbitrary directions of propagation to find that the lower- hybrid drift 
anomalous resistance can be very large indeed, up to a factor of 10 10 or more larger than Spitzer-Coulomb resistance. 

d. Runaway effects. Since the collision frequency is an inverse function of the particle current drift velocity it 
allows for the interesting effect that particles of sufficiently high speed cannot be captured by the electric field. They 
instead get another push and escape as so-called runaway electrons. This effect is known since long and applies to some 
of the anomalous collision processes as well as to Spitzer collisions. The reason is that the wave level saturates yielding 
a constant collision frequency for every mode in question. Hence, fast particles do simply not interact but escape 
like in free ballistic flight. Hence there will always be some particles which behave like freely streaming particles. 
These, when flowing along the magnetic field, constitute a moderately energetic particle beam and may provide a 
seed particle population for further acceleration even in the presence of anomalous collisions. 

e. Other effects. Several other effects have not been mentioned here in relation to anomalous effects. These 
are wave decay processes, generation of radiation, and resonant wave absorption processes in inhomogencitics. It is 
interesting to note that, historically, [3S] in following calculations of [33] proposed that such wave decay processes 
would contribute substantially to anomalous dissipation in subcritical shocks by enhancing the number of waves 
present and thus enhancing the probability of particle scattering by waves. 



Wave decay processes. The first of these belong to the class of wave-wave interaction which are the pendant 
of collisions between particles in the wave picture without invoking particles. Such interactions must satisfy the wave 
momentum and wave energy conservation laws which can be written simply as 

£> a =0, ^> Q (k) = (110) 

a a 

where the index a marks the particular wave mode, and each pair [w a (k),k Q ] satisfies its own linear (or nonlinear) 
dispersion relation. According to the number of waves involved, the smallest possible number is 3, there are three-wave 
interactions, four-wave interactions, . . . with for weak interactions with decaying probability. For the three-wave case 
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one important case is that a strong wave (one that has become very strong due to linear wave growth) such a three 
wave process can lead to wave decay into two other weaker modes of different dispersion due to the relations 

u o — » t0\ + ui 2 , k +ki + k 3 = 0, w >K,W2) (111) 

We do not go into the details of these process; they are rather involved because of the complexity of the three dispersion 
relations which have to be taken into account and which can be different when waves from one branch jump over into 
another branch feeding a wave mode there. These processes have, however, a number of consequences of which four 
are important: 

• They contribute to the excitation of wave modes in a plasma which are not directly driven by an instability 
but result from the decay of an instability-driven large amplitude intense wave, in this way contributing to the 
production of a broad spectrum of turbulence that consists of many different and even possibly weakly damped 
modes in the plasma which otherwise would not exist., when only the decay is stronger than the natural damping 
of the wave. 

• Decay processes limit the intensity of an instability-driven mode and reduce it substantially to the advantage 
of other modes. In this way they weaken the contribution of the particular mode to anomalous collisions while 
they might enhance collisionality by producing a broad turbulent background fluctuation spectrum. 

• By generating other waves in the plasma they may provide a background from which other instabilities may 
grow which are driven by sources which otherwise would not overcome the instability threshold. 

• The broad spectrum produced in plasma wave decay processes may move upstream of the shock and modify the 
upstream conditions in a way not foreseen in the Rankinc-Hugoniot relations. Hence such processes cannot be 
handled in simple plasma modelling of shock wave generation. They can only be investigated by particularly 
tailored numerical simulation techniques. 



Radiation. Only radiation generated from plasma waves is of interest in heliospheric shocks because the densities 
are generally far too low for reaching a substantial emission measure in synchrotron or x-ray radiation. Radiation 
can then be generated only by mode conversion or mode coupling. The difference between the two is that in mode 
conversion an intense high frequency plasma wave propagating up a density gradient gradually transforms into a free 
space mode. 

The more interesting case is a special case of the wave decay process, it is in fact its inversion, when two plasma 
waves join in interaction to inject their energies into a high frequency radio wave that is able to escape from the 
plasma. Radiation production is thus always a process of energy loss that leads to cooling of the plasma. However, 
in the heliosphere this cooling is weak and can be safely neglected, this kind of radiation is in fact degraded to an 
energetically completely unimportant process that has only indicative power. 

The free space modes can be either an ordinary or an extraordinary wave, both propagating above their upper 
cut-off frequencies which for the left-hand polarised ordinary (LO) mode with dispersion lo 2 — uj 2 e + k 2 c 2 is the 
plasma frequency, and for the right-hand extraordinary (RX) mode is a cut-off frequency slightly higher than the 
upper-hybrid frequency to 2 h = 



2 



The following wave-wave processes are of interest in generating radiation: 

• the interaction of two counter-streaming electron plasma (Langmuir) waves into the transverse (T) electromag- 
netic wave, following the symbolic process L + L — > T, where the symbols L, T mean the tuples (wlj^l) an( i 



(wt, /ct), respectively. The energy and momentum conservation equations have been given above in Eq. (Ill |. 
This process produces a transverse wave with nearly perpendicular propagation k^ -C &l in the RX-mode and 
of frequency lo t — 2oj pe , which can clearly propagate above the cut-off in weakly magnetized plasma, 

• the process L + U — > T of interaction of a Langmuir wave with another Langmuir wave that has been scattered 
off thermal ions (i) according to the process L + i — > L' + i* , where the prime indicates the scattered Langmuir 
wave, and the star on the ion the excited ion. The wave interaction in this case causes a lower frequency 
transverse wave still above the plasma frequency but closer to the cut-off of the RX-mode, 

• a process similar to the one under the first item in which harmonic Langmuir waves have been generated in L-L 
interaction. This yields weak higher plasma harmonic radiation at frequencies to < lto pe with harmonic numbers 
1 = 3,... in the RS-T mode with intensity that decreases steeply with increasing harmonic number I, 

• merging of a Langmuir and an ion acoustic wave (S) according to L + S — > T, which produces a transverse wave 
with frequency above but very close to Lo pe . Whether this wave can escape from the plasma depends on whether 
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its frequency is above the cut-off of the RX-mode, which is possible only for high frequency S waves, or whether 
it can tour into the OL mode for which generally the conditions of being selected in the merging process are 
worse than for the RX mode, 

• merging of another electrostatic plasma wave like, for instance, an electron acoustic mode (EA) with a Langmuir 
mode according to L + EA — > T. Since the EA wave is near \uJ pei this produces a T wave at frequency u> ~ § Lu pe 
well above the cut-off still even though still counted as fundamental mode radiation, 

• then, the lower-hybrid mode can also merge with the Langmuir mode. The product of this process can, however, 
not escape in the RX mode as it usually does not exceed the cut-off. On the other hand, it could excite the LO 
mode, 

• finally several nonlinear processes exciting transverse free space electromagnetic waves have been proposed like 
Langmuir wave collapse. Collapse is very attractive because Langmuir waves become trapped in this process 
inside deep density depressions where they bounce back and forth between the walls. Thus in quite a natural way 
collapse generates counter-streaming Langmuir waves in highly localized places, these waves are particularly well 
configurated for merging and escaping in the RX mode. During collapse moreover the internal plasma frequency 
decreases rapidly, causing a decrease of the radiation frequency as well. This allows only the higher frequency part 
of the radiation with frequency just above the plasma frequency of the environment to escape. Radiation occurs 
as intense fundamental but highly bursty radiation at the plasma frequency because at the end of the collapse 
the intensity of the waves explodes and therefore the emitted power also explodes. However, experimentally 
this process could never been proven. It has been replaced by another mechanism known as 'stochastic growth' 
which itself is doubtful as well but very popular. It takes into account the stochastic modulation of Langmuir 
growth in a medium of spatially fluctuating density, i.e. containing a broad spectrum of weak ion acoustic 
turbulence. This causes the growth rate to experience spatial modulations leading to exponential modulation of 
the Langmuir amplitude. Contribution to the intensity comes only from localized places, and thus the volume 
contributing to radiation is a fraction only of the total volume while locally radiation may be rather intense. 



C. Shock particle reflection 

The process of particle reflection from a shock wave is one of the most important processes in the entire physics of 
collisionless shocks, as we have noted already in several places. The mechanism of particle reflection has not yet been 
fully illuminated, however. 

Particle reflection is required in supercritical shocks as, to our knowledge, it is the only process that can compensate 
for the incapability of dissipative processes inside the shock ramp to digest the fast inflow of momentum and energy 
into the shock. Shock particle reflection is not dissipative by itself even though in a fluid picture which deals with 
moments of the distribution function it can be interpreted as kind of an ion viscosity |35j , i.e. it generates an anomalous 
viscosity coefficient which appears as a factor in front of the second derivative of the ion velocity in the ionic equation 
of motion. As such it also appears in the ion heat-transport equation. The kinematic ion viscosity can be expressed 
as 

Mvis = m.A^Anfp ~ P t /2u a (112) 

through the ion pressure Pi and the ion-cyclotron frequency oj C i when replacing the mean free path through the ion 
gyro-radius.. 

In this sense shock particle reflection constitutes by itself a very efficient non-dissipative dissipation mechanism. 
However, its direct dissipative action is to produce real dissipation as far as possible upstream of the shock in order 
to dissipate as much energy of motion as remains to be in excess after formation of a shock ramp, dissipation inside 
the ramp, and reflection of ion back upstream. The shock does this by inhibiting a substantial fraction of inflow ions 
to pass across the shock from upstream into the downstream region. It is sending these ions back into the upstream 
region where they cause a violently unstable upstream ion beam-plasma configuration which excites a large amplitude 
turbulent wave spectrum that scatters the uninformed plasma inflow, heats it and retards it down to the Mach number 
range that can be digested by the shock. In this way the collisionless shock generates a shock transition region that 
extends far upstream with the shock ramp degrading to the role of playing a subshock at the location where the 
ultimate decrease of the Mach number from upstream to downstream takes place. 

Shock reflection has another important effect on the shock as the momentum transfer from the reflected particle 
component to the shock retards the shock in the region of reflection thereby decreasing the effective Mach number of 
the shock. 
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FIG. 17: The two cases of shock reflection. Left: Reflection from a potential well <E>(x). Particles of energy higher than the potential 

energy can pass while lower energy particles become reflected. Right: Reflection from the perpendicular shock region at a curved 
shock wave as the result of magnetic field compression. Particles move toward the shock like in a magnetic mirror bottle, experience the 
repelling mirror force and for large initial pitch angles are reflected back upstream. 



The outcome of the previous paragraphs is that shock particle reflection is of incomparable importance in shock 
formation and in the understanding of collisionless shock physics. On the other hand it is barely understood and can, 
in principle, be treated only by numerical simulations. Before, in the next chapters, coming to discuss those problems 
in greater depth we will present below a few attempts to understand shock reflection. 



Specular reflection of ions from a shock front is the simplest case to be imagined. It requires that the ions experience 
the shock ramp as an impenetrable wall. This can be the case only when the shock itself contains a positive reflecting 
electric potential which builds up in front of the approaching ion. Generation of this electric potential is not clarified 
yet. In the soliton picture the shock potential is related to the density gradient, however, dissipation processes are also 
involved. Understanding its formation requires understanding the entire collisionless shock physics. In a very naive 
approach one assumes that in flowing magnetized plasma a potential wall is created as the consequence of charge 
separation between electrons and ions in penetrating the shock ramp. It occurs over a scale typically of the spatial 
difference between an ion and an electron gyro-radius, because in the ideal case the electrons, when running into the 
shock ramp, are held temporarily back in the steep magnetic field gradient over this distance while the ions feel the 
magnetic gradient only over a scale longer than their gyro-radius and thus penetrate deeper into the shock transition. 

a. Reflection from shock potential. Due to this simplistic picture the shock ramp should contain a steep increase 
in the electric potential A<p which will reflect any ion which has less kinetic energy miV^/2 < eAtfi. This condition 
contains the perpendicular ion velocity component along the shock normal. Since the ion gyrates it depends on the 
instantaneous angle the ion velocity has with respect to the magnetic field at the location of the ramp. In this reflection 
the ion velocity component —Vn — > +Vn simply changes sign. For the gyrating particles this component adds up of 
the normal components of the bulk flow velocity Vjy and the microscopic particle speed v — (v± cos a, v± sin a, vu ), 
with a the gyration angle, projected on the direction of the shock normal n. This yields 



1. Specular reflection 



v N 



VL + vn cos Bn + v± cos a sin b 



'I! 



and the condition for specular reflection can be written as 




+ vii COS&B71 + v± cosasinQsn) 2 < 2eA<j>/m. 
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FIG. 18: Expected particle distributions in the shock foot of a supercritical perpendicular shock which reflects ions back upstream. 
Top: Distribution functions in flow direction. The ion and electron distributions flow with velocity V\ in shock direction. For the 
electrons with their high thermal speed the flow velocity is practically negligible. Reflected ions in the foot occur at velocity — Vi on the 
scale of about one ion gyro-radius. This yields a two-ion beam configuration which is electrostatically unstable. Moreover, ion acoustic 
waves are excited to the right of the reflected ion distribution by the ion-acoustic instability in the positive slope of the hot electron 
distribution. Bottom: Same in the direction along the shock front. Inflowing ions and electrons have only thermal velocities in this 
direction. The reflected ions are accelerated to about twice the inflow velocity in the tangential inflow electric field |E| = |Vi X Bi|. This 
causes an unstable ion beam-ion configuration and a two-stream configuration between accelerated reflected ions and electrons. 



This is a condition on the gyration angle a restricting the gyration phases of the reflected particles. For a gyrotropic 
upstream distribution one can average over all gyration angles from to 7r/2 since only upstream directed velocity 
components reduce the inflow velocity to values below the reflection threshold, obtaining 

(u|| cosOsn - Vl) 2 + -(«|| cos6 B n - v£)v± sin Q Bn + -vj_ sin 2 e> Bn < 2eA< ^ 

This condition must be used on one of the velocity components v\\ , wj_ when determining the number of specularly 
reflected particles from the upstream ion distribution function F™ p (v\\,v±). 

b. Mirror reflection. Another simple possibility for particle reflection from a shock ramp in magnetized plasma 
is mirror reflection. An ion approaching the shock has components v.^ and Vj| = V\ cos(Vi • Bi) along the magnetic 
field. Assume a curved shock like Earth's bow shock. Close to its perpendicular part where the upstream magnetic 
field becomes tangential to the shock the particles approaching the shock with the stream and moving along the 
magnetic field with their parallel velocities experience a mirror magnetic field configuration that results from the 
converging magnetic field lines near the perpendicular point. Conservation of the magnetic moment fi = Tij_/B 
implies that the particles become heated adiabatically in the increasing field; they also experience a reflecting mirror 
force — f^VnB which tries to keep them away from entering the shock along the magnetic field. Particles will mirror at 
the perpendicular shock point and return upstream when their pitch angle becomes 90° at this location [this theory 
has been developed in detail by 33, 55J. We will return to this mechanism in Chapter 6 when describing shock particle 
acceleration. 
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A rough estimate of the marginal upstream pitch angle for mirror reflection can be given from conservation of \x. 
Since the increase in field strength is according to the magnetic gradient across the shock ramp one has roughly at the 
perpendicular shock point B — Bi + (V n B) A where A is the shock width and V n the field gradient across the shock. 
Hence, to lowest order, B = B\ + SB with SB the magnetic compression. This yields for the upstream pitch angle at 
reflection sin 2 oci > (1 + SB/Bx)- 1 . With compression factor SBjB\ ~ 3 particles of upstream pitch angles ol\ > 60° 
will become mirror reflected from the perpendicular shock area due to the action of the mirror force, a condition which 
has to be used upstream in the inflowing distribution if one wants estimating the fraction of reflected particles. This 
requires knowledge of the compression factor, however. The compression factor and the number of reflected particles 
are not independent. Hence, a selfconsistent determination requires developing the full shock theory. This can be 
done only by numerical simulations. 

Of course, the above estimate is very crude. It demonstrates, however, that a fraction of upstream particles can, 
in principle, become reflected from a curved shock surface by mirroring in the converging magnetic field geometry 
around the perpendicular area of the shock. For fast flows reflection will always be located on the nose inflow side 
of the shock. This holds for ions as well as for electrons. Reflection of both sorts of particles has continuously been 
observed at the bow shock as is shown schematically, for instance, in Figure ?? which has been drawn for condition 
in front of a perpendicular shock at a distance inside the foot, roughly within < 1 r C i from the shock ramp. 

c. Consequences of shock reflection. How far the reflected ions return upstream depends on the direction of the 
magnetic field with respect to the shock, i.e. on the shock normal angle For perpendicular shocks the reflected 

ions only pass just one gyro-radius back upstream. Seeing the convection electric field \E y \ — \V?B\ they become 
accelerated along the shock forming a current, the velocity of which in any case exceeds the inflow velocity (which is 
zero in the perpendicular direction) and for sufficiently cold ions also the ion acoustic velocity ci a in which case the 
ion-beam plasma instability will be excited in the shock foot region where the ion current flows. This may generate 
and anomalous collision frequency in the shock foot region. Moreover, since the excited waves accelerate electrons 
along the magnetic field other secondary instabilities can arise as well. 

In quasi-perpendicular and oblique shocks the ions can escape along the magnetic field. In this case an ion- 
beam/ion beam situation arises between the upstream beam and the plasma inflow with the consequence of excitation 
of a variety of instabilities, electromagnetic and electrostatic. In addition, however, an ion beam-electron beam two 
stream situation is caused between the upstream ins and the inflow electrons which because of the large upstream 
electron temperatures probably excites mainly ion-acoustic modes but can also lead to Buneman two-stream mode 
excitation. These modes contribute to turbulence in the upstream foreshock region creating a weakly dissipative state 
in the foreshock where the plasma inflow becomes informed about the presence of the shock. The electromagnetic low 
frequency instabilities on the other hand, which are excited in this region, will grow to large amplitude, form localized 
structures and after being convected by the main flow towards the shock ramp interact with the shock and modify 
the shock profile or even contribute to shock formation and shock regeneration. 

Specular reflection from shocks is the extreme case of shock particle reflection. Other mechanisms like turbulent 
reflection are, however, not well elaborated and must in any case be investigated with the help of numerical simulations. 
In the following subsection we will in passing encounter one such mechanism in phase space hole formation in a two- 
stream unstable collisionless plasma. 

Figure[l8] shows the expected particle distributions for shock reflection of ions in a perpendicular shock right in 
the foot region. The upper part of this figure is perpendicular to the shock in flow direction. The lower part is 
along (tangential to) the shock surface. The different configurations of the distributions in this region can lead to the 
excitation of ion acoustic and two-stream instabilities. Ion beam-ion beam interaction is expected in the direction 
perpendicular to the shock in addition to ion-acoustic instability between the reflected ions and the main electron 
component in the region of positive gradient on the electron distribution. In the direction parallel to the shock 
surface, on the other hand, one expects and ion beam-ion interaction and a two stream instability of the main electron 
component with the accelerated along the shock surface reflected ion beam. Inside the ramp conditions are more 
involved and will be describe in more detail in the context with observations in the respective chapters for both 
quasi-perpendicular and quasi-parallel shocks. 

We note in addition that there is a peculiarity concerning quasi-parallel supercritical shocks. Due to the presence 
of an intense reflected and transformed ion component in the foreshock of a quasi-parallel shock there is a broad 
spectrum of large amplitude low frequency electromagnetic waves which are convected by the inflow stream towards 
the shock ramp, steepen and interact with the shock. These waves are predominantly transverse having components 
tangential to the shock on scales of the order and shorter than the ion gyro-radius. Consequently, quasi-parallel shocks 
remain to be quasi-parallel for ions, in particular for the more energetic accelerated ions, however, for the electrons 
all quasi-parallel supercritical shocks become quasi-perpendicular in the vicinity of the shock ramp transition such 
that for electrons no quasi-paallel supercritical shocks exist. This has the consequence that electrons will become 
reflected and accelerated all-over the shock front independent on its quasi-parallel character. We will later prove this 
statement by referring to observations and simulations. 
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FIG. 19: Left: Electron and ion phase space plots at three simulation times of the interaction of a marginally two-stream unstable 
plasma with a shock ramp. The ramp has been modelled as a density dip (potential wall for electrons) in the centre of the simulation 
box. The box has one space and one velocity coordinate. The fast electron beam current is injected into the quiet ion plasma and [after 
141] . Right: Time history of the elecric field during the interaction with the shock ramp showing the evolution of electron and ion phase 

space holes and their interaction [from [32]. 



2. Hole formation 

As the last item in this section we consider the stability of an high velocity V& > v e electron beam or electron current. 
In the fluid picture we have found that such high speed current flowing through the plasma along the magnetic field 
or (in the presence of steep density gradients) perpendicular to both, the gradient and the magnetic field, will excite 
the fast growing Buneman two-stream instability. It has, however, been predicted early [TBI EEj that currents of this 
strength will undergo a kinetic instability which structures the electron and ion phase space into so-called phase space 
holes which are regions of lacking particles localized in phase space while in real space represent localized electric fields, 
trapped particles and particle acceleration. Such holes have meanwhile been found to exist all-over in collisionless 
space plasmas in relation to spatially localised strong current flow as in reconnection, auroral phenomena, and also 
in shocks [TJ. Since strong currents are expected in shocks as well in the ramp as in the foot, as we have discussed 
above, it is not unreasonable to assume that phase space holes might form under shock conditions as well. 

Hole formation follows a nonlinear interaction known as Bernstein-Green-Kruskal (BGK) mode formation [cf., e.g., 
[6j and is based on the splitting of the phase space distribution function into two components, particles that are 
energetic enough to surpass the electric potential of a localized electric field inside a soliton, for instance, and particles 
of lesser energy that become either trapped or reflected from these potentials depending on the sign of particle and 
potential. Even if the potential is repulsive and ejects, say, electrons from the region, some trapped electrons will 
remain there performing oscillations and become heated up to a certain energy that leaves them still trapped. These 
electrons are in disordered motion and are assuming a high temperature, while the rejection and expulsion of other 
particles from the potential site causes their acceleration. This mechanism is quite complicated and has been treated 
analytically only up to a certain approximate degree in the above cited papers. In order to investigate it one better 
performs numerical simulations. 
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Figure 19 shows the example of a simulation of hole formation in interaction of a marginally two-stream stable current 
(of electron and ion of opposite bulk velocity and same initial temperature T e = Ti with a localized inhomogeneity 
[from |4"T1 |4"2] . The localized plasma inhomogeneity has been modelled as a simple density dip oc — cos 2 (x — Xq) at 
the centre. This is not a shock, it is, however, a potential wall which should reflect one sort of particles, in this case 
electrons. 

The figure shows on the left phase space plots at three different simulation times for electrons and ions when 
due to electron reflection at the potential ramp a strong two-stream instability evolves. The holes formed on the 
electron distribution appear early in the simulation as egg-like distortions of the distribution in the reflected electron 
component. Widening of the distribution indicates the strong heating of the electrons. The holes contain dilute 
trapped electron, and some part of the beam becomes accelerated. At later times the heating becomes violent with a 
strong broadening of the electron distribution when the ion hole starts forming in the lower panels. Strong acceleration 
of a narrow and thus very cool electron beam is also observed in the final state. In addition the holes move along the 
beam, while the ion hole moves in the opposite direction. 

Most interesting is the time history shown on the right in the figure. It shows the initial evolution of many small 
amplitude electron holes moving at fast speed to the right away from the potential ramp. At later times the ramp 
steepens, and the electron holes start interacting with the ion hole which moves slowly to the left. The holes intensify 
and finally can break through the potential ramp to escape to the right where a whole fabric of interfering holes 
evolves. 

The importance of this observation is that two-stream instability can form as a cause of reflection at a potential 
ramp. This is expected for shocks as well. Moreover, the instability causes electron and ion phase space holes to 
evolve and leads to completely collisionlcss heating due to electron trapping inside the holes, i.e. it causes irreversible 
heating and entropy which is needed for shocks, and it also generates a very cool electron beam to escape from the 
holes by continuous acceleration and collimation of cool but fast electrons. This is a very interesting and important 
mechanism which is capable of injecting a fast seed particle population into shock acceleration. 



VI. BRIEFING ON NUMERICAL SIMULATION TECHNIQUES 



The modern age of physics is to a large degree determined by the availability of high speed and high capacity computer 
systems. The use of these computing facilities for performing numerical experiments on collisionless plasmas covers 
now almost half a century of experience. It started with the introduction of Fermi's Monte-Carlo method and 
blossomed after the formulation of the Fermi-Ulam 1961 numerical model approach to cosmic ray acceleration [cf., 
e.g., [34] which was based on nonlinear particle motions in electromagnetic fields. The modern state of the art in 
application to plasma physics has been formulated in several textbooks [e.g., |3] and review articles [e.g., [TOl ITT] . 
Many problems in plasma and in particular space plasma physics with their enormous complexity could not have been 
solved or even attacked without computers and numerical simulations. Also, most of the discussion on shock in the 
following chapters will be based on such numerical simulations which must accompany observation and experiment 
in order to understand what is going on in the shock environment. A brief discussion about numerical methods is 
therefore not only unavoidable but even necessary. 



A. Basic Equations 



Computers are capable to deal simultaneously with the dynamics of many particles as we have already described in 
Chapter 1. The most fundamental approach in numerical simulations is hence based on the full Newton-Maxwell set 
of microscopic equations of which we write down here only the Newtonian subset 

dv. is e s dx 

— — = — (E + VfcXB), — =v 113 

at m s at 

The connection to Maxwell's equations is done by the microscopic definition equations of the space charge p(x, t) and 
conduction current j(x,i) according to 

p(x,t) = ^e s c5(x-x is ), j(x,i) = ^e s v s £(x-x is ) (114) 



In fact, these are the most general equations of a classical plasma consisting of i point-like space charges e s of 
species s with mass m s and momentum m s v is located at time t at location x is (t). The point-like character is 
taken care of by the 5-functions. This whole set is the set of Liouville's equations [which could also formally be 
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combined into one single exact equation in phase space (x, v) by introducing an exact phase space distribution function 
.F(x, v) = ^2 is <5(x — Xi S )<5(v — Vis). Such an equation is know as the Liouville equation]. 

These equations can be simplified depending on the nature of the problem. For instance, when electromagnetic 
effects arc not of interest, the magnetic field will drop out and one uses only Poisson's equation for the electrostatic 
potential and the electric field in Ampere's law. When the plasma is collisionlcss one can use another simplification, 
i.e. replace the Liouville distribution with the one-particle distribution and switch to the Vlasov equation. This 
then produces Vlasov-codes. When one is interested only in low frequency responses of the plasma, the electrons can 
be treated as Boltzmannian electrons, and the electron equations are replaced by the Boltzmann dependence of the 
electron density on the electric potential with the ions being treated as single particles. The corresponding codes are 
the hybrid codes. An even stronger simplification is the assumption of quasi- neutrality when the Poisson equation is 
replaced by the condition N e = N and the electrons are merely an instantaneously reacting neutralising background. 

B. General methods 

Either of the resulting set of equations must be represented in digitalized form in order to be prepared for treatment 
on a computer. One represents the spatial coordinates as a discrete grid and advances the equations over discrete time 
intervals. The choice of space and time steps is prescribed by the necessary accuracy of the result and by the stability 
of the code. The particles in the code are, however, no more point particles but of finite size, i.e. the Delta-functions 
assume bell- function shapes. This has the consequence that particle experience only small angle collisions well suited 
for studying collisionless plasmas. The integration of the equations is then performed in a series on a large number of 
such discrete steps of the finite sized particles the Lagrangean positions of which are used to deposit the charges and 
currents onto the fixed discrete grid points, followed by solving the field equations on these grid points of the many 
spatial cells. This produces self-consistent fields which are used in the next time step and which are smeared out over 
the volume by interpolation in order to construct a field continuum in which the particle orbits are advanced further. 
One is thus working with space and time differences and an interpolation from charge to grid and subsequently of the 
fields from grid to particle. In view of application to the hcliosphcric shocks in the following we briefly discuss the 
discretization only for the electromagnetic case. 

1. Units 

Numerical methods require that all quantities including coordinates space and time are pure numbers. One thus needs 
to normalize them by introducing constant reference values, a density No, charge-to-mass ratio e/mo, velocity of light c 
(for instance, other choices are Alfvcn velocity and so on). Time is then normalized, for instance, to plasma frequency 
expressed in these reference units u~ 1 = (e 2 N a / e m )^ ^ , space in inertial lengths c/u> p , the electric potential for 
instance in m c 2 /e, the magnetic vector potential in m c/e. This choice of units is one of many possibilities only. 

2. Discretization 

Both the field and particle equations must be discretized. The idea of discretization is quite simple. One returns in 
history to the time just one step before Newton and Leibniz. Differential quotients become quotients of differences, 
higher order differential quotients become the corresponding powers of quotients of differences, mixed differential 
quotients become products of quotients of differences. The only trick is to assign the results to some location inside 
the difference interval, not necessarily the centre (!), and to do this properly Also, time runs only in one direction: 
forward. Applying such a scheme one arrives at recursive equations which can be solved on a sufficiently powerful 
computer. 

Let us assume we reduce the electromagnetic set of equations to the electromagnetic wave equation for the vector 
potential component A. This must be written in difference form 

Jv 2 (A<!) + 2A<*) + A<-*>) - ^2 (1 + /?£> 2 V 2 )(A(§) - 2A<*> + A<-*>) - 

The superscripts indicate the time levels of the various terms, the ad hoc parameter is introduced to modify the 
dispersion at short wavelengths, and V 2 is defined as 

a ^j+i,m 2Aj m + Aj—i tTn Aj m +i 2Aj m + Aj m —\ I 1 i 

~~ ~Kx 2 Ay 2 ' D 2 ~ Ax 2 Ay 2 
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The transverse current density in the Coulomb gauge is jt = j — Vd(j>/dt. since V • Jt = 0, the additional equation 
V 2 {d(j)/ dt) — V • j must be solved. The electrostatic potential (j>^ is to be taken a full time step (1), while j^) is at 
half time step (^). The former is obtained from V 2 ^ 1 * 1 = —p^\ Finally, the fields follow from 

E« = -V<£« - Ail) ~ Aa \ B (D = 1 V x (Ad) + A<$>) 
Y At 2 y ' 

Fields, charge and currents are defined at the centre of the cells. These equations must be solved with a so-called 
Poisson solver. In addition one needs to specify appropriate boundary conditions at the boundaries of the simulation 
box. 

In a similar way one discretizes the equation of motion of the particles. For this one defines h — eAt/m and obtains 
a centred form of the equation of motion as a recursion equation with unknown v*- 1 ) 

v(i)=v(5) + / l (E( 1 ) + v' 1 ) x B«) 

In order to determine B = |(v^i) -f-v*-^) the former equation is solved implicitly taking the scalar and vector products 
of the former equation with B*- 1 )). This yields up to terms second order in At the expression 

,(§) = (\ - ^^i*fF.iv(i) - / '~ 



v^> = ^1 - y B J + h ( E + v xB) + y|ExB + BB • v ( 3- 

for use in the expression for v' 1 ) only. The set of equation obtained is then ready for computing. 

However, there are two ways of computing, so-called explicit and implicit techniques. In the explicit technique one 
solves the equations as they are obtained after discretization. In such an approach the internal errors will necessarily 
grow and at some stage become unstable such that the calculation must be truncated. One can artificially introduce 
some kind of damping for these growing error modes in order to suppress them. Justification for this is discussed in the 
literature. The implicit technique solves the problem by calculating backward in time |17j such that the error modes 
decay away when time runs negative. We do not describe this technique here. It suffices to note that in this approach 
the dangerous unstable short-wavelength error modes disappear by definition and become partially eliminated from 
the system. The most efficient ways of calculation are the combinations of both methods. 

All these methods work on a fixed grid. In application to shock ramp research the grid mash has to be chosen 
refined enough for resolving shock structure. Recently [25] developed a self-adaptive technique which takes care of 
the steepening and narrowing of a shock front in order to resolve its substructure. This is a significant progress in 
shock simulation technique. 

We do not go into detail of the various methods and refinements of simulation techniques. Those readers who 
are interested and prepare by themselves for doing simulation research in collisionless shocks we rather refer to the 
mentioned basic literature on the techniques of numerical simulation. 
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